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INTRODUCTION 


Under the conditions typically prevalent in highly loaded transonic fan 
compressor stages, the linearized, small-perturbation approximations to the 
equations of motion cannot be expected to be descriptive of the flow in the 
vicinity of the blades. Thus, recourse is made to the numerical solution of 
the complete nonlinear system of equations, as discussed in Reference (l) in 
connection with the b 1 ade- to-b 1 ade program. However, sufficiently far from 
the blade rows, the amplitude of the flow disturbances will decay to acoustic 
levels and the linearized, sma 1 1 -perturbat ion approximations will be descrip- 
tive of the far-field. Therefore, an intermediate region in which both analy- 
ses are valid should exist at some distance from the blades. The inlet and 
discharge stations of the b 1 ade-tO“b 1 ade computational domain can serve as the 
interfaces betv/cen the near-field (numerical) and far-field (acoustic) analy- 
ses. The present far-field analysis is formulated with respect to an infinite 
duct model, namely, all outgoing waves should propagate without reflection. It 
differs, however, from conventional inlet duct analyses in that the signal may 
begin with an arbitrary transient, associated with the deviation of the assumed 
initial data in the near-field from the periodic solution which is sought as the 
asymptotic limit in time. Therefore, the acoustic analysis must recognize that 
a transient signal will occur during startup and that a simple harmonic time 
dependence, which is the usual basis of inlet duct acoustics, cannot be assumed. 
The analysis should allow the transient to radiate outward without reflection, 
and should be capable of identifying the attainment of a periodic solution by 
the growth of discrete harmonic components in the solution. 

INTERFACE WITH NEAR-FIELD SOLUTION 

The procedures for defining inlet and discharge stations in terms of axial 
boundaries of the computational domains shown in Figure (l) have been described 
in References (l) and (2). The present analysis pertains to the flow upstream 
of the inlet station and downstream of the discharge station, where the passage 
is assumed to be an annular segment of an infinite cylindrical duct. With re- 
spect to the analysis of Reference (l), the streamsheet is assumed to have con- 
stant radius, r, and thickness, b. However, the circumferential boundaries of 
the passage are not necessarily restricted to a single blade-to-blade passage 
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FIGURE I. BLADE-TO-BLADE COORDINATE SYSTEM AND GRID NETWORK 



(as in Reference l), but should include the angular interval 2 tt/AN, where 
AN = (In the special case AN=0, which is usually avoided in practice, 

use of a single blade- to-bl ade passage is again permissible.) Use of the pre- 
sent analysis with the computer code described in Reference (2) will, therefore, 
requi re addi t tonal storage and retrieval of data along the interfaces between 
the numerical near-field solution and the present acoustic far-field solution, 
in the same manner as is performed along the interface between domains h and 5 
of the near-field solution (see Reference 2). 

Derivation of the boundary conditions on the near-field (finite difference) 
solution, at arbitrary stations referred to as the inlet and discharge boundaries 
of the near field, was presented in Volume 1 (Reference 1). The acoustic analysis 
is formulated in the absolute frame of reference, and accordingly the boundary 
point analysis from Reference (1) will be restated here in absolute coordinates. 


As discussed in Reference (l), the boundary point solution is obtained from 
a reference plane me t hod-of-character i s t i cs procedure.* A pair of (approximately) 
two-dimensional wave motion characteristics and a stream surface characteristic 
are identified at each grid point on the inlet and discharge boundaries, as 
shown in Figure (2). The compatibility relations, 


^ ± 
dt 


pa 


dV 

n 

dt 


= - YP ( 


r36 




''e V dm^ 


apply on the wave motion characteristics, 

dm rdO 

v~ 

m 0 

which are depicted by tfe lines AO and CO In Figure (2). The equation, 


S = constant 

applies on the stream path characteristic, 

dm _ rdO _ 


*As pointed out in Reference (1), translation of the reference plane effec- 
tively transforms a rotating frame of reference back to an absolute frame. Thus 
the numerical results are independent of the coordinate system in which the 
equations are stated. 
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which is depicted by the line BO in Figure (2). In addition, if the flow is 
Isentropic 

C/p = constant 

also applies on the stream path characteristic, BO. Finally, the circumferen- 
tial momentum equation is written as: 


!!§ 

dt 


8H 

rae 


+ T 


as 

rae 


- ? 


( 5 ) 


The inlet flow is assumed to be isentropic and i r rotat iona 1 ; thus: 


S. = 0, or p./pT = p /p^ 

I I I -oo ^ -oo 

and 

= 0 


giving: 


( — 
'at ' 


Inti educing the small 
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disturbance approximations, 


( 6 ) 

(7) 

( 8 ) 


P == P ■ P-cx> P. 


V = 


V - V << V 
mm m 




« V 
0 m 


etc . 


the compc: 1 1 b i 1 i ty equation on the outgoing wave, CO, at the inlet becomes: 


(9) 

( 10 ) 

( 11 ) 
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J ^-oo -oo m, r'c 
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-oo rn 


1 av^ V ^ ^ 

* f drbx X 

" 2 ^YP-co W . ■'y Pc + ?b ^ 
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( 12 ) 
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where the subscript i denotes the inlet value at a new time, t+At, i.e,, at 
point 0 in Figure (2), and subscript c denotes the value at point C. Recall 
that r and b are assumed constant for rn ^ m.. Since Vq can be evaluated at 
all points 0 along the inlet, from integration of Equation (8), and all 
variables at point C can be obtained by Interpol at ion of known (current) 
data in the near^field, the right-hand side can be regarded as a known function 
of (m. , 6 , t) , i.e.: 


P; 


p., 


a 

— oo 


\1' 

m. 

I 



im. , 


e, t) 


(13) 


Q. (n.^0, t) taken as the boundary condition to be placed on the up- 

stream far-field solution. As will be shown in the following section, this 

condition is sufficient to determine p. , and thus p.. can then be obtained 

I I nij 

from Equation (13), and p, from Equation (6). As indicated above, \/q is ob- 

1 D • 

tained from Equation (8), completing the inlet solution at the new time t+At. 


Similar procecfures at the discharge boundary lead to the relations: 


Pd^pd = Pfa/pb 




'at >— 


i*L 

rao 


+ V c), 

rae m ^'d 


(11*) 

(15) 


Pd Pco ^co \ = Poo 4 % (mj, e, t) (16) 

d 

Equation (l6) represents one boundary condition on the downstream far-field. 
Integration of Equation (15) provides a second boundary condition on the down- 
stream far-field, which, is required since it cannot be assumed to be i r rota t i ona 1 . 


Ve Rd ("’d' (’7) 

d 

Due to the finite work input by the rotor, the downstream reference con- 
ditions (p^, a^) are also not necessarily the same as the upstream reference con- 
ditions this case, the far-field analysis provides the value of 

p^ and Equations (lM, (l6) and (17) complete the discharge station solution at 
the new time tHAt. 
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With these preliminaries in hand, attention is now focussed on the acous- 
tic far-field analysis. In the following discussion the primes on the perturba- 
tion variables will be dropped as the acoustic solution is understood to per- 
tain onMy to perturbation quantities- Furthermore, the subscript ( 
will be used to denote the reference conditions for either boundary, that is, 

X ± and all variables without subscript will refer to perturbations with 
respect to the reference state, for example, p = p - p^ and V = V - V^. 

FAR-FIELD ANALYSIS 

The problem of interest is the acoustic radiation from the considered blade 
rows under steady-state operating conditions. However, development of the nu- 
merical solution for the near-field includes generation of a transient flow field 
associated with the deviation of the assumed initial conditions from a periodic 
solution, which is sought as an asymptotic limit in time. Therefore, analysis 
of the far-field must necessarily recognize that a transient signal will occur 
during startup and that a simple harmonic time dependence, which is the usual 
basis of engine acoustic analyses, cannot be assumed to be descriptive until 
the asymptotic limit is approached. The analysis should allow the transient to 
radiate outward, without reflection, and should be capable of identifying the 
attainmenl of a periodic solution by the growth of discrete harmonic components. 

Use will be made of a harmonic type analysis to develop an acoustic formula- 
tion of the pressure field in which frequency effects are superimposed and the 
time history obtained by integration. This procedure leads to a convolution 
integral (i.e., a type of Duhamel integral) in the time domain, which is more 
directly compatible with the numerical data in the near-field than a result in 
the frequency domain. (The near-field data consists of values of the pressure, 
flow velocity, etc., at a discrete series of grid points and at fixed time in- 
terva Is.) 

Since the present effort is addressed toward a cascade formulation, the 
governing equations are written in a two-dimensional, Cartesian coordinate sys- 
tem. It is noted that the two-dimensional problem could be described by a so- 
lution of the wave equation alone, were it not for the fact that the time de- 
pendent force distribution on the blades is capable of producing a convected 
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vorticity field (as well as a corresponding entropy field, which, however, is 
not relevant to the present problem). It will be seen that the convected vor- 
ticity field does not contribute to the acoustic pressure field, per se, as dis- 
tinguished from the irrotational component of the velocity field which is di- 
rectly coupled with the acoustic pressure. Therefore, the velocity perturba- 
tion field downstream of the discharge boundary and upstream of the inlet 
boundary (see Figure l) is characterized by the sum of an irrotational (Vx=0) 
velocity vector and a solenoidal (V*=0) velocity vector which satisfy 
the linearized conservation equations. Note that the inlet flow has been as- 
sumed to be irrotational, therefore, in the upstream far-field; however, 

the analysis will be developed with respect to the more general case pertaining 
to the downstream far-field. 

The linearized conservation equations for the perturbation variables can 
be stated as* 


where : 



•V = 0 

(18) 


0 

(19) 

dp 2 dp _ 

dt dt 

^0 dS 
C dt 

V 

(20) 



(21) 

11 

< 

3 

+ 


(22) 

V. = I V + 

2 1 TI 2 


( 23 ) 

7 X = 0 


{2k) 

7 X ^2 = ^ 

0 

(25) 


*These represent linearized versions of Equations (13)» (l^)» (15) and (17) of 
Reference (l), with r = constant and b = constant. 
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= 0 

(26) 


U = M a 

( 27 ) 


0 00 


and 

A+ u 

dt dt 0 9m 

(28) 


Since the streamsheet has been assumed to be cylindrical in the far-field, 
with b << r, the conventional Cartesian coordinates (x,y) will be used in place 
of the meridional coordinates (m,0) : 

(29a) 
(29b) 

Substituting Equation (26) into (l8) gives: 


X “ m 
y “ re 


dp 

dt 


+ 


2 

a 

o 



= 0 


(30) 


Substituting Equations (24) and (26) into the curl and divergence of Equa 
tion ( 19) , respect ively, gives: 


4z + — v^p = 0 (31) 

dt 1 Po 

and : 

^ (VXV 2 ) = 0 ( 32 ) 


Equations (I 8 ) and (19) or ( 30 ) and ( 31 ) yield the wave equation for the pres- 
sure: 


1 

2 

a 




0- 


o 

Therefore, solutions for p are purely radiative. 
( 31 ) can be also combined to yield the wave equation, 
purely radiative. However, Equation (32) infers that 
convect i ve . Therefore : 


(33) 


Since Equations 
solutions for 
solutions for 


( 30 ) and 
are also 
are purely 



dt 


= 0 


(3M 
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and: 


^ + J_ 7p = 0 
dt 


Solutions to Equations (33), (3^) and (35) can be expressed in terms of a 
Fourier integral representation: 


p(x,y,t) 




V (x,y,t) 


> = Ze 

f n 


[v (x,y,t)J 


- 1 a y 


( , .\ 
p a A (o)) 
o o n 


hi.) 

^c^i.) 


) a e 
f o 


2tt 


(35 


(36 


' / s'' 

V (x.y.t) 




2 


Vf. (x,y,t) 
"2 



D (to) 
n 


E„(») 


i(o(t-x/U ) , 

i o dto 

/ 2ir 


(37 


Consider first the downstream boundary on which (see Equations 16 and 17) 

.2n 


<id(y>t) = (Pd Po ®o V^/(po ^o) 


Cl 


(38 

(39 


Thus, the boundary data can be expressed as: 




Q(y,t) 

tR(y.t)j 


)= 


“lay 
n' 


CO 

f 




Q»(») 


V i tot dto 
) ® 2 tt 




where ( ) denotes the time to frequency transform, as above, while ( )^ indi- 

cates transformation from spatial location, y, to spatial harmonic, n, which 
will be discussed later (cf. Equation 55). (The subscript d has been dropped 
since the ensuing development can apply to either the inlet or discharge 
boundary. ) 


(AO 
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Substitution of the above integral forms into Equation (2^), (26), (32), 
(3^), (38) and (39) gives: 


C 

n 


E 

n 


U a 
o n 


D 

n 


B 

n 


a 6 
o n 

0)“U Br, 

o n 


+ B + D 
n n 


n 



C 

n 


+ E 

n 


n 


(41) 

(42) 

(43) 

(44) 

(45) 


The value of A , and thus B , C , D and E , can thereby be expressed In terms 
n n n n n ' 

of the transform of the boundary data, viz.; 


w-U B 
o n 


o)^+(a -U„)6, 


+ U a a 
o o o 




U a R°] 
on n 


(^6) 


The propagation coefficient is determined by substituting the integral rela- 
tion for pressure (Equation 36) into the wave equation (Equation 33): 




-w M ± (o) 
o 


no o 


o o 


(^7) 


where the + sign refers to downstream propagation (from the discharge boundary) 
and the - sign to upstream propagation (from the inlet boundary). 


Equations (46) can now.be substituted into the corresponding integral rela- 
tion (Equation 36) to obtain, after some manipulation, the pressure perturbation 
on the discharge boundary in terms of the specified boundary data. 


P a 
^o o 


pd(y.t) = 2 

0-m:) 


2 e 
n 


lay 

n' 


a a 
n Q 


CO 

/ 


(03-Sj ( 


n' 'a a n 
n O 


4. M ^ 

* % V ^ 27 


(^8) 

11 


I 



where 


S = (( 1 )^ - (1-M^))^ 

n no o' 


(A 


Consider now the inlet boundary where V 2 = 0. The boundary condition here is 
stated as (see Equation 13): 


Q. (y.t) = (Pj - Pq v^_)/(pQ 


(5' 


KJ 

Since = 0, D = E =0, and R becomes part of the solution, i.e., it can- 
not be specified as a boundary condition, consistent with the near field anal- 
ysis described in the previous section. Thus, is only a function of 0 ., viz: 


-o 

\ 0) - 


is 


where 


Q| (y,t) = 


r 

n In 


(u) 


^ io3t 

2tt 


( 5 - 


Again, with some manipulation, an equation for the pressure perturbation on 
the inlet boundary is obtained in terms of the specified data; 


P{(y,t) = 


P 

^O O 


(1-M^) (1+M^) 


00 

-i« y 
E e " 

n 2 2 / 

a a / 
n o 


(w-S^) (o) 


( 5 : 


\ n® IWt dw 
+ M ) Q e 7 ^ 
o n' n 2 tt 


The selected representation of the solutions, Equations (3^) and (37), and 
boundary conditions, Equations (^O) and (52), as Fourier series in the y direc- 
tion is appropriate for enforcement of the periodicity boundary condition per- 
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taining to spatial variations in this direction. The coefficient is de- 
fined accordingly: 


“n = 


liin 


nAN 


( 54 ) 


where Y is the fundamental wave length of the stage (double cascade) configura- 
tion (i.e, Y = 2iTr/AN). In addition, the fact that the boundary data is speci- 
fied at a discrete number of grid points, say N, on the boundaries implies that 
the Fourier series can only include N terms; n=0, 1, 2...N-1. Since the 
distance y to each point can be written as mY/N, where m=0, 1, 2...N-1 also, 
the Fourier series can be expressed in the standard Discrete Fourier Transform 
(DFT) notation: 


f ) 

r 


Q(y.t) 

S / = S 

1 

R(y.t) ^ 

1 

^ / 

1 



N-1 

Z 

n=0 


Q“(t) 

R°(t) 


2iri nm 

~~1T— 


(55) 


The invcije DFT is then: 


/ ^ 

(f{t) 

n 

k 

1 i 

/ \ 
i 

R»(t) 

\ j 

m=0 

R„(t) 

^ m / 


27t 1 nn/N 


(56) 


For cor.ipii ta t i ona 1 purposes, the following non-dimensional parameters are 
evident in Equations (48) and (53): 


Q = (jj/aa 


T = aa t 
o 


(57) 

(58) 


where 


a = 2 tt/Y = a /n 
n 


(59) 
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Then, on the discharge boundary: 


where: 



D a 2 / 

p”U) 

PQ^O 1 

■ (1-M2)n2 J 

T 

n 

- (1-M^)n 

O 

Q°(fi) 

-o 

= aa Q (o)) 

o n 

R°(fi) 

n 

-o 

= aa Q (w) 

o n 

inlet 

boundary : 


(fi Q°(n) + n R° (n)e'^'^) || 




( 60 ) 


( 61 ) 

( 62 ) 

(63) 


P°(x) 


O O 


(1-M^) (1+M )n^ 




(64) 


Then, on either boundary: 

N-1 

p (t) = T. e 

^ n 

n=0 


•2'irinm/N o 


Pn(x) 


(65) 


The Integrals of Equations (60) and (6^) can be transformed back to the time 
domain by the following convolution. On the discharge boundary: 


,v Q^(t)^+ J^(t) R^(t)) 

n o o n n n n 


( 66 ) 


where * denotes convolution in time. For n > 1: 


H°(t) = - 

o 


CO 

/ 


_ X \^T dQ 

n (fi“T ) e — 

n 2"^ 


(67) 
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J°(x) 


1 


n(l-M^) 

o 


oo 

/ 


/ ^ N dQ 

' 5T 


and, far n=0: 


H°(x) 


6(t) 


J°(x) = 0 




where 6 Is the Dirac delta function. 


Note that for t ^ 0 with n 0: 


H°(x) j 6(x) 


p °(0 


•t p a^Q°(x) 

2 o on 


On the Inlet boundarcy: 


P°(x) 




where, for .n > 1 


K°(x) = 1 

n 2/, „2 


n"(l-M") (1+M^) 

o o 


oo 

/ 


(J2-T ) (n+M T ) e'^'^ — 

n on 2 tt 


and, for n=0: 


( 68 ) 

(69) 

(70) 
(71) 


(72) 

(73) 

( 71 *) 

(75) 


15 


( 76 ) 


K°(t) 


j5(t) 




while for t 0 with n 0: 


k®(T) = j6(t) 


P^(t) = “ p 

2 o o^n 


( 77 ) 


(78) 

(79) 


The integral in Equation (67) is evaluated by dividing the non-dimensional 
frequency range into three regions, viz: -» to to and to ». 

In the middle region standard DPT techniques are applied, while in the outer 
regions, the integral is accurately approximated in terms of an ex- 

ponential integral of the form 


+0O 



— o 


for which closed form expressions are available. The integral in Equation 
(68) can be expressed in terms of a Bessel function: 


oo 

f (n- §. AW (8,) 

— OO 


2 4 

where c denotes n(l-M^) and A(t) is the Heaviside step function. The evalua- 
tion of the integral in Equation (75) is carried out in the same fashion as 
that in Equation ( 67 ). 
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The desired solution for the pressure perturbation on the interface with 
near-field is finally accomplished by observing that Equation ( 65 ) can now 
be expressed in the form: 


m 


2 

o'^o 


N-1 


n=0 


-2TTinm/N < 




K°(t) QO(t) 


R°(t) 


) ( 82 ) 


on the discharge and inlet boundaries respectively. Use of Equation (55) for 

P (t) , Q„(t) and R_(t) and similar DFT expansions for H , J and K then leads 
m ' ^ m m m 

to the convolution: 


p a" N-1 


* Vn") * ” Vn<‘> 




K (t) » Q (t) 
n m-n 


(83) 


on the discharge and inlet boundaries respectively. Thus, a double convolu- 
tion over both time and distance (in the y direction) is required. The func- 
tions* snd R^(t) , therefore, represent the point sources of time-varying 

strength v;hich are aligned along the considered boundaries, with spatial reso- 
lutions consistent with the number of grid points specified. 

The functions H^(t) , d^j^(t) and *^j^(t) are the duct response functions, 
defined by: 


*Bear in mind that Q-^.(t) denotes either Q. (y ,t) at the inlet station or 
Qd 3t the discharge station, while denotes R^ (y^,t). 
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> 


\ 

H^(t) 

m 

N-1 

H°(t) 

m 

n=0 

.°(t) > 

K (t) 
m 


, ^n<'> 

J 


J 


■ 


^6(t) 




N-1 
+ Z 
n=1 


-27Tinm/N 


H“(t) 




K°(t) 


“27rFnm/N 


( 84 ) 


Finally, it should be noted that only the real part of p (t) is of interest- 

m 

Equations (67), (68) and (75) infer that H^(t) and K*^(x) are purely real 
o n n ' 

functions and J^(t) is a pure imaginary function. The boundary data Qj^(t) 
3nd fife obviously all real numbers. Therefore, it is concluded that 


J (t) = Z - i J°(t) sin 
m n n N 


27rnm 


I 0/^\ • 27rnm 

Z L (t) s m —r^ — 

n n N 


(85) 


where 


L°(t) 

n 


1 


n(l-M^) 

o 


oo 


( 86 ) 


I.e. , the solution can be carried out entirely in terms of real functions, 

NUMERICAL EXAMPLE 


A group of subroutines including an efficient Fast Fourier Transform 
routine has been developed, for use in conjunction with the blade-to-blade 
computer program B2DATL, described in Reference (2), to carry out the indica- 
ted convolutions required to solve Equation (83) numerically. Results have 
been thus far limited to test cases with a simple harmonic input signal. For 
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example, in one of the calculations the discharge station was assumed to be 
divided into 8 intervals covering a total circumferential distance of 0.1 m. 

The selected reference (average) Mach number of the discharge flow was 0.8. 

The flpw was assumed to be irrotational so that only one input function Q(y,t) 
was required, and the second function R(y,t) could be considered as a response 
function (that is, it was calculated from (i(y,t)). The input function 
Q(y,t) = cos (fix - Zirn/N) , with Q(y,t) = 0 for x < 0, was selected for this 

case, where x = 2TTNa^t/Y, (n = 0, 1, 2, N-1) N = 8, Y = 0.1 m., 

a^ = 10^ rn/sec, and>i2 == 1. The input function Q and response function R at 
n = 0 are plotted in Figure (3). It should be noted that the response function 
is initially out of phase with the input function because of the assumption that 
Q = 0 for X < 0. However, the effect of the transient at x = 0 dies quickly, 
and after about 1/3 millisecond the response function closely approximates the 
input function and Indicates the desired harmonic solution Is being approached 
asymptotically. The non-dimensional perturbation pressure is plotted in Figure 
{h) , The complete history is shown for the point n = 0, whereas the history 
of the points n = 1 and 2 is only shown at early times, where a difference in 
amplitude as well as phase exists. The values of p^ and p^ are not plotted at 
later times (that is, after about 1/3 millisecond) since the pressure solutions 
at the various grid points only differ noticeably In the phase angle corres- 
ponding to the input function and the differences in amplitude asymptot i a 1 1 y 
decay. However, the complete numerical results for n = 0, 1, 2 7 are in- 

cluded in Appendix B as the first test case. 

DESCRIPTION OF COMPUTER CODE 
Overview 

Program RAFFT (Real Acoustic Far-Field Theory) , 1 1 sted in Appendix A, is 
designed to compute the real part of the acoustic far-field of a pair of inter- 
acting cascades in an infinite duct. This program has not as yet been coupled 
to program B2DATL which computes the near field of the cascades. It currently 
consists of two parts. The first part, which is made up of the main program 
RAFFT, the subroutine GETQR (which calculates the input functions Q and R) and 
the function routine TR00T, is at present used solely for the test purposes; 

It must be replaced by a calling routine in B2DATL which calculates the input 
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PRESSURE PERTURBATION 




functions Q and R from the boundary data In accord with Equations (38), (39) 
and ( 50 ). The remainder of the program consists of seven subroutines and func- 
tion routines: AJ1DX, C0NV0L, FFT, FFT2, GREENS, P0LY and TC0NV, Only sub- 

routine CjCNV 0L Is called from the main program; all other subroutines and func- 
tion routines are called from C0NV0L. Subroutine C0NV0L is called by the state- 
ment: 

CALL C0NV0L (P, Q, R, NY, IT, DTIME, MAXC0N, CMARK, RH0O, C, Y, FIRST, 
NDIMY, NDIMT, PU) 

which requires all variables as input parameters except P and PU which are the 
only results. P(IY, IT) is the discharge station perturbation pressure, as a 
function of discrete spatial location space, lY and time, IT, and PU(lY,IT) is 
the corresponding inlet pressure. 


Input Data 

The input parameters Q(IY,IT), R(IY,IT), NY, IT, DTIME, MAXC0N, CMARK, 
RH0O, CY, Y, FIRST, NDIMY and NDIMT, required for the evaluation of P(lY,IT) 


through 

the statement CALL C0NV0L (, 

. .), are d.^fined as follows; 

NAME 

DIMENSION 

DEFINITION 

Q 

(8,1001) 

Near-field discharge or inlet bound- 
ary condition, (P^lPo^^V^) /p^a^. 

See Equations ( 38 ) and (50). 

R 

(8,1001) 

C i rcumferent i a 1 perturbat i on velocity, 
V^/a^. See Equation (39). 

NY 


Number of grid points in y-direction. 

IT 


Cq.rrent index. 

DTIME 


Time step; seconds. 

MAXC0N 


Maximum number ofvpoints used for con- 
volution in time (< 1001, with present 
dimension of Q and R) . 

CMARK 


Reference Mach number, M^. 

RH0O 


Reference density, p^; kg/m^. 

C 


Reference speed of sound, a^; m/sec. 
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NAME 


DIMENSION 


DEFINITION 


Y 

0MEGAO 

FIRST 

NDIMY 

NDIMT 

Subroutine and Fun 

The following seven subroutines a 
portion of the program necessary to de 


Fundamental wavelength of the cascade 
(27rr/AN) ; m. See Equation (5^). 

Input frequency (non-dimensional), 

See Equation (57) • 

A logical parameter which must be set 
to **FIRST = ,TRUE,'* before subroutine 
C0NV0L is called for the f irst time. 
This avoids repetition of calculations 
that need to be done only once. Sub- 
routine C0NV0L automatically sets 
"FI RST=. FALSE. ‘‘ after the initial call. 

Dimension of the arrays Q and R in the 
Y (circumferential) direction, presently 
8. (NY <_ NDIMY is necessary.) 

Dimension of the arrays Q and R in 
time, presently 1001. 

tion Sub-Programs 

d function routines encompass the main 
ermine the acoustic far-field solution: 


NAME 

C0NV0L 


AJIDX 


DE SCRIPTION 

ViMien called the first time, C0NV0L calculates a number of 
parameters used through the rest of the program, calls sub- 
routine GREENS to calculate the Green*s functions H and 

m 

Thereafter, it calls subroutine TC0NV to perform the 
convolutions, and finally it evaluates Equations ( 83 ): 

NY 

P(IYJT) =*^— (i: [Q(IY,IT)*H(KY,IT)+R(IY,IT)*J(KY,IT)]) 

KY=1 

where * denotes convolution in time. 

Evaluates the Bessel function J1 (X)/X using a 7th order 
polynomial approximation. 
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NAME DESCRIPTION 

FFT Converts a complex array into two real arrays, and then 

calls FFT2. 


FFT2 


A fast Fourier transform subroutine written in assembly 
lanquage*. 


GREENS 


PC^LY 


Calculates the Green's functions H , J and K In accord 

m m m 

with Equation (84) for m = 0, 1, 2. ...NY-1. In this 
routine the expression f (W) =W=v^^ - 1 occurs where W = Q 
(Equation (57)). The sign of the square root is chosen in 
such a way that f(W) is symmetric in the imaginary part 
and antisymmetric in the real part, with f(o)=+i. To cal- 
culate (t) we evaluate the integral in a finite interval 
from -WCUT to +WCUT using the fast Fourier transform, where 
the cut-off frequency is calculated as WCUT=Y/C/DT IME . The 
contribution from the tails outside the interval (-WCUT ,+V^CUT) 
is just 1/2 the delta function, plus terms of second order 
(i.e., containing the factor 1/WCUT**2). 

Evaluates a polynomial, using as few multiplications as 
possible. 


TCCInV Forms the convolution of two arrays A and B, calculating 

only one new value 

•^max 

S. = E A(j)*'B(i-J) (for one index i=IT) 

' j=i 

where i =m?n(i, MAXC0N) ^ and MAXC0N is the maximum number 
•'max 

of points used for convolution. 


*0nly suitable for use on CDC 6600 computer systems. -However, similar FFT 

routines are generally available for other systems. 
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Flow Chart of the Program 


The interconnection between the subprograms are represented in the 
schematic below: 

Part 1 (Generates input functions and prints output. Should be replaced 
by a calling routine from B2DATL.) 



FFT2 


P0LY 







Test Cases 


The Input function (J (t) = cos (fi t - 2tt •^)*A(t), where t = ZirNa t/Y 

n N o 

is non-dimenslonal time, and N=NY is the number of points in the y-direction, 
was used for a test example. In this case the response function, 

rather than input function. 


To check R^(x), it was compared with its asymptotic value. Two cases 
must be distinguished. 


If n/l-M^ 

5 




^n, imag^^^ ~ 

1 

[n - T (fi)] sin 
n 

Q T 


n(1-M) 

u 

n^m 

'^n, real “ 

1 

[fi - T^(fi)] cos 

^ T 

6 

n-m 

n(l-M) 

with “ 

1 

CM 

(1-M^) 



If < n/l-M^, 

substitution of T^.^ = - 

i|T 

nl yields: 


1 

s i n T + T 

n 

COS 

^ «n-m 

n(l-M) 


1 

[Q cos Q T ~ T 

n 

s i n 

n t] 6 

n-m 

n ( 1 -M) 


( 87 ) 

( 88 ) 

( 89 ) 

(90) 

(91) 


Two test examples were run, both with m=1 , and with Q = 1,0 in the first 

case and Q = 0,5 in the second case. The first test case corresponds to the 

numerical example discussed in Section IV, Due to the factor 6 , R = 0 

n-m’ n 

except for n=l. 


In the first example Q 
= 0.5 < = 0.6. 


K* > n/l-M =0.6. In the second example 


As shown in Appendix B, R^ approaches its asymptotic value reasonably 
well in both cases. 
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The print-out in Appendix B includes: 

• a list of the values of all parameters used 

• a comparison of R(x) with its asymptotic value 

• the input functions Q(lY,IT) and R(IY,IT) 

• the Green's functions H^(t) and J^(x) and 

• the resulting pressure P(IY,IT). 



APPENDIX A 


PROGRAM LISTING 


PROG^?AM RAFFl (INPUT, OUT PUT , TAPERS) 

DIMENSION P(6, 1001 ) #0(6, 1001 1001 ) ,PU(8, icon 

LDGICAL FIRST 
dimension ASYNPKS) 

complex EPOuER, TRUOT , t2PI#Rl »R7,Ewl ,Ert7 
COMMON/NSTtP/NSUP 
EPOWtR(X)=CEXP(C*^PLX(0, ,X) ) 

* 

* READ AND PRINT PARAMETERS 

10 READ 11 0,nnME,NY,NT ,MAXCON,CMARK,RHOO,C» Y^wZERD 
no FORMAHFIO.O, 51 10,«F 10, 0) 

IF(NT.EO,0) STOP 

PRINU20,DT1ME,NY,NT,MAXCON,CMARK,RHOO,C, Y,InZERO 
120 FORMAT (*i'pARAmET£RS USED*/ 

** DTIME »*fcl0,3/ 

** NY **110/ 

** NT **110/ 

** MAXCON «*I10/ 

*« CMARK «*E10,5/ 

** RHOO «*E10,3/ 

** C **E10.3/ 

** y **£10, 3/ 

** QMEGAO **tl0.3/ ) 

ft 

NSTEPC5 
NDIMYsS 
NDIMTbI 001 
FIRST*, TRUE* 

ft 

11 CAUL SECOND(TO) 

PRINT IM 

lia F'JRMATC///* TEST RISTAR APPROACHES ITS ASYMPTOTIC VALUE*/) 

DO 12 1T*1#NT 

ft GET NEVi values OF Q AND R AT TIME T FOR NY DISCRETE VALUES OF Y 

call G£ T0R(0, Ri it ,DlIMt ,NY,CMARK,C# Y,NDIHY,MDIMT,>nZ£RO) 

12 CONTINUE 

ft 

PRINT 115 

115 FORMATC*! time input R(71ME) •) 

DO 15 IT*1 ,NT,NSTEP 
TlMEaCIT*l)*DTIME 

15 PRINT l25,TlMf , (R(IY,1T), lY=lrNY) 

125 FORMAU 1X,E 10.3,FR.3,7F7,3) 

PRINT 116 

116 FORMATC*! TIME INPUT O(TIME)*) 

DO 16 lT*l,NT,NSTtP 

UME*( I T-1 )*DTIME 

16 PRINT 12S, time, (Q(1Y,IT)#1Y*1,NY) 

CALL SECONOCTI) 

T JSED=Tl-TO 
PRINT 310,TU3E0 

ft 

DO 20 1T*1,NT 

* calculate P(Y,T) for discrete values of V 

* BY convolution #.1TH GREFNS FUNCTIONS 

CALL CONVOL (Pr D#R# NY, IT , DT IMF, MAXCON, CMARK, RHOO, C, Y#F I R5T# 

ft N01MY,NDJMT , PU) 

ft 

ft PRINT SELECTED VAUTS OF THE RESULTING PRESSURE P 
1F(M00(IT-1,NSTFP).NE ,0) GO TO 20 
1F(IT,E0,1) PFINT 130 

130 F0RMAT(*1 time pressure P(TIMU for DISCRETE VALUES OF T*) 

IFCIT.EQ.l) CALL SECONO(Tl) 

T1ME*( IT-1)*DTIME 

PRINT 1«0,TIv.E, (P( IY, IT), IY*1,NY) 
lao FORMATC IX, <?E10. 3) 

20 CONTINUE 

CALL SEC0NDCT2) 

TUSED*T2*T1 
PRINT 310,TUSED 

310 FORMATC//* TIME USED **F7,3* SECONDS*//) 

PRINT 150 

150 FORMATC*! TIME UPSTREAM COMPONENT OF PRESSURF*) 

DO 50 1T=1,NT,NSTEP 
TIME*CIT-1)*DTI^E 

50 PRINT 1U0,TIME,(PUCIY, n ),1Y»1,NY) 

END 
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SUBROUl INE GETO^ (Q,R,IT,OTU^E,NY,CMARk,C, y,fJDIMYfN01MT,rtZERO) 


DIMtNSlUN 0 (NDIMY,^DImT),H(N0IMY,N0IMT) 
heal Ji(ioon,j7nooi) 

complex tgiSTAKClOOUf U7STAP(100J)fRlSTAR,R7STAR,RR(8) 

complex ASYPP7, trout 

COMMON/NSTfcP/NSTEP 

PI»3,1«IS9P6535898 

T«0PI»2.*PI 

Ts(lT-l)*DTlME 

AtPMACaTrtClPI*C/Y 

WO»ALPriAC*wZeRO 

WOTs^'OaT 

A«S0RT(1,-CPARK**2) 

OD 10 lYslrNY 

0(IY, IT)=CUSCWOT-TttOPl*(IY-l)/NY) 

10 CONTINUE 

Jl (IT)=(l i+CMAPK)*AJlDX(ALPHAC^A*T) 

OlSTARnn = .5*CEXPCCMPLX(0.,ft0T)) 

R1STAR=0, 

DO 20 JT=l,n 

R1STAR = R1STAR + J1 (JT)*CJlSTAR(lT-JTfl) *(0,, 1.) /PI 

20 continue 

IF(HOaUT-l,N5TEP).NF,0) GO TO ?S 

ASYMPTs ( 1 ./( 1 .-C^'ARK) ) * ( rt ZERO- TROOT ( 1 , CMA RK , rt ZERO) ) 

A * CEXPCCPPLX(0.,^.OT)) 

PRINT 125rlT,M|STAR,ASYMPT 

125 FORMATC* tire STEP*l«* PlSTAR **2F7,3* ASYMPTOTE s*2F7,5) 
25 CONTINUE 

DO 30 1Y*1,NY 

R< IY» IM=REAL( CEXP ( CMPLX ( 0 w-T rtOP I * ( 1 Y- ! ) /N Y ) ) *R 1 ST AR ) 

30 CONTINUE 
RETURN 
END 


complex function TKllOT(lYrCMAf?K,W) 
Sart**2-lY**2* (1,-CMARK**?) 
IfCS.GE.O.) TROOTiSORT (S) 

IF (S.LT.O.) TROtn?(Ow-l.)*SORT(-S) 

RETURN 

END 
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SUBROUTINE CON VOL ( P # Q , R . N Y , I T , DU '^E , MAxCON, C^ARK , RHOO , C f V r F I RST , 

• NUIMYrNDlMT,PU) 

* 

t GIVEN THE INPUT FUNCTIONS tHY,T) = ( P + RHOO* C *U ) / ( RHOO *C ) AND R(Y»T)=V, 

* FIND THE PRESSURE P A'l XsO THROUGH CONVOLUTION IN TIME AND SPACE 

* P(IYrlT) = RHOO*C*(1/N) « SUH(OVKR KY FRO^■ 0 TO N-l) OF 

* ( 0(KY,1T) CONV H(1Y-KY,H) ♦ R(KV,IT) CONV J(IY-»^Y,IT) } 

A 

DIMENSION P(NDIWY,NDI^T ) , (J ( NO I MY , ND I MT ) r R ( ND 1 MY , NDI MT ) 

DIMENSION 0C0NV(S,S),RC0NV(8» fl) 
real HR(1001,S)fJR(1001#5) 

REAL PUC'JOI^^Y,NDIMT ) ,OKCONV(6r8),KR(100l,5) 

logical FIRST 

COMMON/PI /PI, PI half, Tr^OPI 

DATA PI/3, l«15R?hS35698/ 

* 

A IF THIS SUBROUTINE IS CALLED FOR THE FIRST TIME, 

A calculate the greens FUNCTIONS HR AND JR 
A AND THE NON-DIf^ENSIOUAL TIME STEP PT 
1F(, NOT, FIRST) GO TO 10 
PIHALF=PI/2, 

TaOPI=2,*PI 

rhoocn=rhooaC/ny*c 

OT=DriME*C*TivOP‘l/Y 

NY?a,\Y/2 

NY21XNY2-M 

NY22=NY2+2 

f01MF = NDIi *t.OIM/2 

CALL SECOND(TO) 

call GREENS(HR, JR,DT,MAXCON,NY,CMARK,P,NOJMP,KR) 
call SECONO(Tl) 

TUSEDsT 1-TO 
PRINT 310,TUSED 

310 FDRMAT(//a TIMf USED =*F7,3 a SECONDS*//) 

firsts. false. 

A 

A PREPARE A table of TivtE CONVOLUTIONS 
A (SAVES COMPUTER TIME flY USING SYMMETRIES) 

10 DO SO IYX1,NY 

DO 20 KY31 ,NY2l 

call TCONVCUCIY, 1 ) ,HR( 1 ,KY) , IT,DT,NY,MAXC0N,0C0NV(1Y,KY) ) 

CALL TCONV(Q(Il ,1)»KKn,KY),IT,0T,NY,MAxC0N,0KC0NV(IY,KY)) 

20 CONTINUE 

DO 30 KY1CNY22,NY 
KY-NY+2-KY1 

0<C0NV( 1Y,KY1 )=QKCUNV( I Y,KY) 

30 0:ONV( I Y,KY1 )=C3CONV( 1 Y, KY) 

A 

RCONV(lY,l)sO. 

RC0NV(IY,NY21):'0 
00 liO KYs2,NV2 
KY1SNY+2-KY 

CALL TCONV(wn Y, 1 ) , JR( 1 , KY-1 ) , n,DT ,fgY,MAXCON,RCONV(I Y,KY) ) 

^0 RC0NV(IY,KY1 )=.RC0NV(1Y,KY) 

SO CONTINUE 

A 

DO 70 IYS1,NY 
SUMO=0, 

SUMRxO, 

SUMKsO, 

DO bO Ky* 1,NY 
IKYslY-KY+1 

IF(IKY,LT,1) ikyciky+ny 
SUMU aSUMQ+OCONVCKY, IKY) 

SUMR*SUMR + PCGNV(KY, IKY) 

SUMKsSUMK + OKCONV(KY, IKY) 

60 CONTINUE 

PJ( lY, J DsRHOOCN*SUMK 
70 P(IY, IT)=RHOOCN*(SU‘-’Q + SUMR) 

RETURN 

END 



SUBROUTINt GREENS (HR, JR, D1 ,MAXCON, Y , CHARk , P , , KR) 

A generate the greens functions HR(T) and JRCT) 

A fOR 1 « Or OTr 2*DTr (MAXCON-I)*OT 

A And for R s 0, l, NY/2 (for H), R « l, 2, NY/2-1 (FOR J) 

A 

REAL HR(lOOlrS), JR(1001,3),COSJKC6),SINJK(B) 
real KR(l001r5) 

COMPLEX S,P(l) 

COMMON/PI/P I , PI half, T rtOPI 

* 

A P, WHICH IS NOT DEFINED AT THE BEGINNING# IS USED 
A AS DUMMY STURAGE FOR THE FAST FOURIER TRANSFORM 

* 

A FIND MAXIMUM N ^HICH FITS WITHIN SPACE OF P 
Nalb364 

10 IF(N,LE,N01FP) GO TO 20 
N»N/2 
GO TO 10 

* 

20 CM«J ,-CMARKa^» 2 
A«SfJRT(CM) 

WCUT=].^UPI/DT 

DwswCUT/N 

NY2=NY/2 

NY21=NY2+ I 

NYMJaNY-1 

NY2M1SNY2-J 

NasN/4 

A TABLES OF SINES AND COSINES 
DO 25 r=l,NY 
ARCxTwOPI*(l-n/NY 
SINJK(I)=SIN(ARG) 

25 COSJK(1}=COS(ARG) 

A 

A CALCULATE HI 

1000 FORMAT(1X,M0.3#F9.3,6F7,5) 

1160 FORMAT(a1 time GREFNS FUNCTION HR*) 

PRINT 1100 
DO 30 Irl,N 
H« ( I-l ) *Dw 
S03WAA2-CM 

IF(SQ,GE,0.) SxSORT(SQ) 

IF(SQ,LT,0, ) S = CMPLX(0*,-SORT (-SQ)) 

30 P(I)«(rt-S)*w/CM 

A 

CALL FFT(P,N) 

A 

DO 40 1=1 ,NU 

A HN IS EVEN, TAKE T^TtE REAL PARI OF FOURIER TRANSFORM 
P( I)=Rt AL(P( I ) )ADr;/PI 
A DELTA FUNCTION 

IF(I.EU.l) PCI)=,5/0T 
40 CONTINUE 

« 

A COSINE SERIES 

DO 70 I=1#MAXCUN 
DO 60 Jl =1 #NY21 
JxJl-1 
SUM=0, 

A HO * 1/2 A DELTA FUNCTION 
IFd.EO.n SUM*.5/0T 
DO 50 K=1 ,NYMi 

1F(KAI .GT.N4) GO TO 50 

INOEX=MOD( J*h,NY)+l 

1F(I.GT,1) SUm=SUM+CUSJK( INDEX) * K*P{K*1 ) 

IFd.EU.D SUMeSUM + COSJKCINDEX) a P(l) 

50 CONTINUE 
60 MR(I#J1)=SUM 

IF(MUD(Ir50).EQ.2 .AND. I.NE.2) PRINT 1180 
T»(I-I)*DT 

160 PRINT lOOOrT, (HRd, J), J*1,NY21) 

70 CONTINUE 

A 

A JN (FOR NXI) IS THE BESSEL FUNCTION M * J 1 ( A*T ) / ( A* T ) 

A WITH A 3 SURT(1-M**2) 

1190 F0RMAT(*1 TIME GREENS FUNCTION JRa) 

PRINT 1190 
DO 60 I=1,N4 
T*dd)AOT 

80 P(I)=CMARK*AJ1DX(A*T) 
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* » 


SJNE SERIES 

03 190 I=l,«AXC0N 
00 100 J*WNY2MJ 
SUM=0. 

00 90 KaUNYMl 

IF(K»I •GT.N4) GO TO 90 

INDEX»><0D(J*K,NY)*1 

SUM»SUMtSINJK(INDEX) • K*P(K*1 ) 

90 CONTINUE 
100 JR(I* J3=SUM 

ir(MOO(I,50).E0.2 .and, l,NE.2) PRINT 1190 
T«(I-l)*OT 

PRINT lOOO.T, (JR(lrJ)r J»lfNY2Ml) 

190 CONTINUE 

« 

* CALCULATE Kl 

1200 F0RMAT(*1 time greens FUNCTION KR«) 

PRINT 1200 
03 130 I»lrN 
w* ( 1 * 1 ) * 0 h 
Sa*rt**2-CM 

IFCSO.GE.O.) S*SQPT(SCJ) 

IFCSO.LT.O.) S=CMPLX(0.,-SnRT(-SO)) 

130 P(l)a(rt-S)*(K + r;MAPK*S)/(CM*( 1 .♦cmARK)) 
call FFT(P.N) 

DO 1«0 I»l,N« 

P(I) = REAL(P(n)*Of-/PI 
IFCI.EO.l) P(I)=.5/0T 
190 CONTINUE 

00 170 I«1,KAXC0N 
DO 160 Jl=l,NY21 
JxJl-l 
SUMsQ. 

IF(I.EO.l) SUMe,5/DT 

DO 150 K=1,NY^^1 

IF(K*I .GT.N9) GO TO 150 

INDEX=moD(J*K,NY)+1 

IFd.GT, 1 ) SUM = SUM + C0SJK(INDEX)*K*P(K*1 1 

IF(1 .Ea, 1) SUM=SUM+COSJK ( INDEX) *p( I ) 

150 CONTINUE 
160 KR(I,Jl)=SUM 

IFCMOO(I,50),EG,2 .AND. I,NE.2) PRINT 1200 
Tx(l*n*DT 

200 PRINT lOOOrT, (KRCI.J), J=1,NY21) 

170 CONTINUE 

• 

RETURN 

END 



SUBROUTINE ICON V ( OR , HJ # 1 T . DT , NY , MAXCON# RESULT ) 


* 

A CONVOLUTION IN TIME 

* OF Q(IY#IT) rtlTH H(JY,IT) 

A OR R(IY,IT) WITH J(JY,ITI 

* TO calculate unf new point 

* 

DIMENSION ORUIrHJCl) 

resulted, 

MAXTeMAxCON 

IFUT.LT.MAXCON) MAXT*IT 
Je(IT-1)*NY+1 
DO 10 K»1,MAXT 
RESULT=RESULT+QR( J)AHJ(K) 
10 JeJ-NY 

RESULTeRESULTaOT 

RETURN 

END 


FUNCTION AJIDX(X) 

* BESSEL FUNCTION Jl(X)/X 

DIMENSION CH7),C2(7),C3(7T 

data (Cl(I),l = l;7)/,b,-,b6aa‘^Q6S, ,2l095b75r-,039S«269, 

* , 00 '1 93M9,-, 00051 76!,, 0 0 001 109/ 

DATA (C2( 1 ) , 1 = 1 , 7)/,7 97S?^yb6. , 00000156/, 01 6‘3‘?66 7,, 0001 71 05, 

* ,002995 11, ,00113653/, 000? 00 35/ 

DATA IC5(I), 1 = 1 ,7>/-2, J5bi99£i9, ,l?099bl 2, , 0 OOOSbSO , • , 00657 87 9 , 

* ,00079596,, 00079629, -,O00?9l 66/ 

AX«ARSCX) 

IF(AX.GT,3,) GO TO 10 
AJlDXsPOLVCCl, (X/3, )**2,7) 

RETURN 

10 F1*P0LYCC2,3./AX,7) 

THETAiaAX+POLY(C3,5,/AX, 7) 

AJlDXsFl*COS(ThETAl )/SQR1 (AX)/AX 

RETURN 

END 


FUNCTION POLY(A,X,N) 

* evaluate Tr(E POLYNOMIAL A Cl ) ♦ A ( 2 ) * X + * . . ♦ A (N ) * X • * ( N- 1) 
DIMENSION A(l) 

P0LY=0.‘ 

XPOhER=1 , 

DO 10 iBlfN 
POLY=POLY + ACn*xPOWtR 
10 XPOwrR*XPDwER*X 
RETURN 
END 


SUBROUTINE FF T ( XRE AL. , NS I G ) 
real XREAL(I) 

IFCNSIG.LT.O) GO TO 10 

NaNSIG 

NSPACE»2 

CO TO 20 

10 Ns*NSlG 

NSPACEb-2 

20 CALL FFT2(XREM,XHEAL(2) ,N,NSPACE) 
RETURN 
END 




IDENT 

FFT2 


NX5 

xq 


ENTRY 

FFT2 


RX7 

XUXb 

LlOO 

sxo 


83 


RXO 

X0 + X5 


S\Mi 


BO 


NXl 

X7 


SB5 


83-B7 

L30 

SB5 

Bbtpq 


AXO 


1 


SA2 

R1 fBq 


SB5 


80 


SA3 

BUBS 


SB6 


XO 


SAq 

B?4Bq 


sxi 


85 


RX6 

X24X3 


to 


86,H7»FFT23 


SA5 

B240S 

LllO 

SBU 


B3-Bq 


RX2 

X2-X3 


SB5 


B3-B5 


SA6 

A2 


SA2 


Bueq 


RX7 

xq + xs 


SA3 


BUBS 


RX3 

XUX2 


SAq 


B2 + Bq 


Rxq 

xq-X5 


NX 7 


X2 


SA7 

Aq 


SA5 


R2 + B5 


RX5 

X0*xq 


NX6 


X3 


RX2 

X0*X2 


SA7 


A3 


RX6 

X3-XS 


SA6 


A2 


Rxq 

XUxq 


NX7 


Xq 


$Ab 

A3 


NX6 


X5 


RX7 

X24Xq 


SA7 


A5 


Sbq 

BbtPS 


SA6 


AU 


SA7 

AS 


LT 


86, Bq, LI 10 


LT 

Bq,e3,L30 

L120 

SSA 


flq + e? 


SB5 

Bq-B3 


SB5 


P64BS 


bxi 

•XI 


SA2 


BUBq 


SBq 

B6-R5 


SAi 


BUBS 


LT 

BS,bq, L30 


SAq 


B2fBq 


SBq 

pq4R7 


NX7 


X2 


SA2 

SO 


SA5 


B?4B5 


LT 

Bq,BS,L?0 


NX6 


X3 

Lao 

SBU 

BO 


SA7 


A3 


SXS 

Bb 


SA6 


A2 


AX5 

1 


NX7 


xq 


SBb 

XS 


SXO 


B6 

L50 

SB5 

Bb4Bq 


NX6 


XS 


SA2 

eupq 


SA7 


AS 


SA3 

8UHS 


SAb 


Aq 


SAq 

B?4Hq 

L130 

AXO 


1 


RX6 

X24X3 


IXI 


xuxo 


SA5 

B24BS 


Pi 


XI ,U 30 


RX7 

X2-X3 


LXO 


1 


SA6 

A2 


s«q 


Bq4P7 


SA7 

A3 


IXI 


XUXO 


RXb 

Xq + XS 


SB5 


XI 


sea 

B64PS 


GE 


8S,Bq,Ll 10 


RX7 

xq-X5 


LT 


Bq,Bb,Ll20 


SAb 

Aq 


SAl 


AO 


SA7 

AS 

FFT2 





LT 

Bq,P3,LS0 


SAO 


A1 


EQ 

B6, 07, LlOO 


SAl 

AO 



SAl 

A1 


SBl 

XI 



SBU 

B7 


SAl 

AO+1 


RXb 

xi*xi 


S92 

XI 



SAl 

AUl 


SAl 

AO + 2 


FXb 

X64Xb 


S35 

XI 



SA3 

ONE 


SAl 

AO+i 


SAb 

CO 


SB-q 

XI 


L60 

BXO 

XI 


SAa 


Oq 


RXb 

X3-X6 


MX2 


1 


BX7 

XO 


SA5 


L60 


NXl 

Xb 


8A5 


B3 


SA7 

SD 


LX2 


57 


EQ 

L30 


PX7 


X3 

S 

DATA 

9.SP7379909S9773qbf-S 


8X6 


-X2*X5 


DATA 

1 ,91 7q7S97M 0 7033U -a 


PL 


xq,Ll 0 


DATA 

3,P3q9SlP7S71 39SS9K-q 


BX6 


X?4XS 


DATA 

7.b69903}P7q270US3fc-q 


exa 


-xq 


DATA 

1 ,S3 3'/H0 ♦ «b2Pq7bSM-5 

UO 

LX3 


32 


DATA 

3,0b79Sb7b296S97b3E-3 


SAb 


AS 


DATA 

b, 1 3sepqbq9isqq7sqf. -i 


NXO 


HS, X3 


DATA 

1 ,?271S382HS7199?(>.f-2 


Px2 


xq 


DATA 

?,qSqi2?eS229l22MHt-2 


S87 


xq 


DATA 

q.9(ib767q327qib01«E-2 


DX7 


X?*X7 


data 

9,801 71 q03295b0b0?E-2 


SAl 


854S 


DATA 

1 .9S0903?201bl2M27F-i 


SB 3 


X7 


DATA 

3,62bS3q323bSC0«i77t-l 


SB6 


X7 


DATA 

0.707l0fa7Bll865q7S 


EO 


Lqo 

ONE 

DATA 

1.0 

L20 

SA3 


CD 

CO 




RXA 


X2^X1 

SO 




RX7 


xa*xo 


END 



RXS 


X3< XO 





RX6 


X3*XI 





Rxq 


xq-x5 





RX6 


X64X7 





APPENDIX B 


PRINTOUT OF TEST CASES 


PARAMETERS 

USED 

OTIME 

s 

.asoE-os 

NY 

M 

6 

NT 

U 

201 

KAXCON 

m 

101 

CMARK 

u 

.BOOEtOO 

RHOO 

s 

.lOQF^Ol 

C 

s 

.lOOEfO^i 

V 

s 

.lOOEfOO 

OMEGAO 

« 

.100£>»01 


TEST RJSTAR approaches ITS ASYMPTOTIC VALUE 


TIME 

STEP 

1 

RISTAR 

• 

0.000 

.143 

ASYMPTOTE 

s 

1,000 

0,000 

TIME 

STEP 

b 

RISTAR 

B 

-,31b 

.757 

ASY'^PTOTE 

a. 

, 707 

.707 

TIME 

STEP 

11 

RISTAK 

X 

-.963 

.928 

ASYMPTOTE 

X 

.000 

1.000 

TIME 

STEP 

lb 

RISTAR 

a 

-1.520 

,520 

ASY'flPTOTE 

3 

-.707 

,707 

TIME 

STEP 

21 

RISTAR 

X 

-1.596 

-.257 

asymptote 

S 

-1,000 

,000 

TIME 

STEP 

2b 

RISTAR 

X 

•1,065 

-.969 

ASYMPTOTE 

= 

-.707 

-.707 

TIME 

STEP 

31 

R.ISTAR 

3 

-.150 

-1,209 

ASYMPTOTE 

= 

-.000 

-1,000 

TIME 

STEP 

3b 

RISTAR 

X 

.701 

-.635 

ASYMPTOTE 

3 

,707 

-.707 

TIME 

STEP 

til 

RISTAR 

3 

1.068 

-.057 

ASYMPTOTE 

= 

1,000 

-.000 

TIME 

STEP 

46 

RISTAR 

B 

,801 

,667 

asymptote 

3 

.707 

,707 

TIME 

STEP 

SI 

RISTAR 

B 

,101 

.983 

ASYPPTOTE 

B 

,000 

1,000 

TIME 

STEP 

5b 

RISTAR 

3 

-,59b 

,660 

asymptote 

= 

-.707 

,707 

TIME 

STEP 

61 

RISTAR 

B 

-.674 

-.023 

asymptote 

3 

-1,000 

.000 

TIME 

STEP 

bb 

RISTAR 

B 

-.560 

-.690 

ASYMPTOTE 

3 

-,707 

-.707 

TIME 

STEP 

71 

RISTAR 

X 

.098 

-.940 

ASYMPTOTE 

3 

-.000 

-I .000 

TIME 

STEP 

7b 

RISTAR 

S 

,736 

-.603 

asymptote 

3 

, 707 

-.707 

TIME 

STEP 

61 

RISTAR 

B 

,944 

.113 

ASYMPTOTE 

3 

1,000 

-.000 

TIME 

STEP 

6b 

RISTAR 

C 

• 577 

,782 

ASYMPTOTE 

3 

,707 

.707 

TIME 

STEP 


RISTAR 

B 

-.159 

1,004 

asymptote 

3 

,000 

l.OOO 

TIME 

STEP 

96 

RISTAR 

C 

-.635 

,636 

ASYMPTOTE 

3 

-.707 

.707 

TIME 

STEP 

101 

RISTAR 

B 

-1,052 

-.110 

ASYMPTOTE 

B 

-1,000 

,000 

TIME 

STEP 

106 

RISTAR 

3 

-,67<? 

-,80b 

asymptote 

3 

-,707 

-.707 

TIME 

STEP 

111 

RISTAR 

= 

.094 

-1.046 

ASYMPTOTE 

3 

-.000 

-1 ,000 

TIME 

STEP 

11b 

RISTAR 

X 

,610 

- , 687 

asymptote 

S 

,70 7 

-.707 

TIME 

STEP 

121 

RISTAR 

3 

1 ,Obb 

.063 

ASYMPTOTE 

X 

1,000 

-.000 

TIME 

STEP 

12b 

RISTAR 

B 

.717 

.772 

asymptote 

X 

,707 

.707 

TIME 

STEP 

131 

RISTAR 

B 

-.031 

1.029 

asymptote 

3 

.000 

1 ,000 

TIME 

STEP 

13b 

RISTAR 

B 

-.743 

,669 

ASYMPTOTE 

B 

-.707 

,707 

TIME 

STEP 

141 

RISTAR 

X 

-l.OlO 

-.045 

ASYMPTOTE 

C 

-1,000 

,000 

TIME 

STEP 

146 

RISTAR 

3 

- ,68 3 

-.7*45 

asymptote 

3 

-.707 

-.707 

TIME 

STEP 

151 

RISTAR 

X 

,039 

-.996 

ASYMPTOTE 

3 

-.000 

-1 ,000 

TIME 

STEP 

156 

RISTAR 

S 

,726 

-.660 

ASYMPTOTE 

S 

,707 

-.707 

TIME 

STEP 

Ibl 

RISTAR 

■ 

,974 

,065 

ASYMPTOTE 

K 

1,000 

-.000 

TIME 

STEP 

16b 

RISTAR 

B 

.637 

.751 

ASYMPTOTE 

X 

,707 

./07 

TIME 

STEP 

171 

RISTAR 

S 

-.064 

.991 

asymptote 

S 

,000 

1,000 

TIME 

STEP 

176 

RISTAR 

X 

-.762 

, bua 

ASYMPTOTE 

X 

-.707 

.707 

TIME 

STEP 

161 

RISTAR 

a 

-.993 

-.068 

ASYPPTOTE 

X 

-1,000 

,000 

TIME 

STEP 

166 

RISTAR 

B 

-.637 

-.775 

asymptote 

X 

-.707 

-.707 

TIME 

STEP 

191 

RISTAR 

a 

.102 

-1,011 

asymptote 

X 

-.000 

-1,000 

TIME 

STEP 

196 

RISTAR 

a 

.791 

-.657 

ASYMPTOTE 

X 

.707 

-.707 

TIME 

STEP 

201 

RISTAR 

a 

1.026 

,065 

ASYMPTOTE 

X 

1.000 

-.000 


time 

fl 

Input 

0.000 

R ( TIME ) 

,101 

.143 

.101 

.,000 

-.101 

-.143 

-.101 

TIl^E 

0. 

input 

1,000 

O ( TIME ) 

,707 

-.000 

-.707 

-1,000 

-.707 

-.000 

• 707 

# l 25 E-oa 

-.316 

.312 

.757 

,759 

,316 

-.312 

-.757 

-.759 

. 125 E -04 

,707 

1,000 

,707 

,000 

-.707 

-1,000 

-,707 

-.000 

,250 t-oa 

-.963 

-.025 

.928 

1.337 

,963 

.025 

-.928 

-1,337 

.250 E -04 

,000 

,707 

1,000 

,707 

-,000 

-.707 

•1,000 

-.707 

,37 b £-0« 

-1.520 

-.707 

,520 

1,442 

1,520 

,707 

-,520 

-1,442 

•375 E .04 

-,707 

,000 

,707 

1,000 

.707 

-.000 

-.707 

•1,000 

, soot.-oa 

-1.596 

-1,310 

-.257 

,94 b 

1,596 

1.310 

.257 

-,94 b 

, S 00 E -04 

-1,000 

-.707 

,000 

,707 

1.000 

,707 

• 000 

-.707 

• 62 bt »0 *l 

- l ,06 b 

•1,439 

-.969 

,066 

1,065 

1,439 

.969 

-.068 

.625 E -04 

-.707 

-1,000 

-.707 

,000 

,707 

1,000 

.707 

,000 

, 7 bO£-OU 

-.150 

-.961 

-1.209 

-.749 

,150 

,961 

1.209 

,749 

•750 E -04 

-.000 

-.707 

-1,000 

-.707 

.00 0 

.707 

1,000 

• 707 

,87 bt -0 a 

.701 

-,09 b 

-.035 

•1,066 

-,701 

.095 

,635 

1,066 

.875 E -04 

,707 

-,000 

-.707 

•1,000 

-.707 

,000 

,707 

1,000 

, lOOE-Oi 

1,060 

.715 

-.057 

-.796 

-1,068 

-.715 

,057 

,796 

• lOOE -03 

1,000 

.707 

-.000 

-.707 

- I.OOO 

-,707 

.0 00 

,707 

• 1 1 3 t »05 

.501 

1,052 

,667 

-.081 

-.801 

-1,052 

-.687 

.061 

• I 13 E -03 

,707 

1,000 

,707 

-.000 

-.707 

•1,000 

•,707 

.000 

, 12 bE -03 

• 101 

.766 

.963 

,623 

-•101 

- , 766 

-,903 

-,623 

•125 E -05 

,000 

.707 

I.OOO 

,707 

-,000 

-.707 

-1,000 

-.707 

• I JOE - 03 

-.596 

,060 

• 680 

,902 

♦ 596 

-.060 

-.680 

-,902 

,130 E - O 5 

-,707 

,000 

.707 

1,000 

,707 

-.000 

-.707 

-1,000 

• 150 E -03 

-,67« 

-.634 

-.023 

,602 

.874 

,634 

,023 

-,602 

•150 E -03 

-1,000 

-.707 

,000 

,707 

1,000 

.707 

-,000 

-.707 

■ 163 E -03 

-.580 

-.904 

-•696 

-.064 

• 580 

,904 

,698 

,064 

. l 63 t -0 i 

-.707 

-1,000 

-.707 

,000 

, 707 

1,000 

,707 

-,000 

• l 7 bt -03 

.098 

-.596 

••940 

.,734 

-•098 

,596 

,940 

• 734 

•175 E -03 

-,000 

-,707 

-1,000 

• ,707 

,000 

,707 

1,000 

,707 

t 1 88 E -03 

.738 

,095 

-.603 

.,948 

-•738 

-.095 

.603 

,946 

•100 E - O 3 

.707 

-.000 

-.707 

-1,000 

-.707 

,000 

• 707 

1,000 

,200 t -03 

.9 U 4 

,740 

.113 

.,586 

-.94 4 

-.740 

-.113 

,588 

•200 E -03 

1.000 

,707 

-,000 

-.707 

-1,000 

•,707 

,000 

.707 

•21 3 E -03 

.577 

,962 

.762 

,145 

-.577 

-.962 

-,762 

-.145 

•213 E -03 

,707 

1,000 

,707 

-.000 

-,707 

•1,000 

-.707 

,000 

• 22 bE - o 3 

-.159 

,598 

1,004 

,622 

.159 

-.598 

-1,004 

-.022 

•225 E -03 

.000 

,707 

I.OOO 

,707 

-,000 

• ,707 

-1,000 

-.707 

• 238 E -0 J 

-,05 b 

-.140 

,638 

1,041 

,835 

• 140 

-,630 

-1,041 

•236 E -03 

-,707 

,000 

.707 

1,000 

• .707 

-.000 

-,707 

•1,000 

• 2 bOE-Oi 

-1.052 

-,022 

-, 1 10 

• 666 

1,052 

.022 

.110 

- • 666 

•250 E -03 

-1,000 

707 

,000 

,707 

1,000 

, 707 

-.000 

-.707 

« 263 E -0 5 

-.672 

-1,046 

-.606 

-.095 

• 672 

1.046 

,006 

.095 

•263 E -03 

-.707 

•1,000 

-.707 

,000 

,707 

1,000 

,707 

• , 000 

,27 bE -03 

.094 

-.673 

-1.046 

-,806 

-.094 

.673 

1,046 

,606 

•275 E -03 

-.000 

-.707 

-1,000 

-.707 

.000 

,707 

1.000 

,707 

• 28 QE >-*0 3 

.810 

.066 

-.687 

•1.058 

-.610 

-.006 

,667 

1,056 

.260 E - O 3 

,707 

-.000 

-.707 

-1,000 

-.707 

,000 

,707 

l.OOQ 

•300 E -03 

1,066 

,796 

• 063 

-,709 

•1,066 

-.790 

-.063 

,709 

•300 E -03 

1,000 

,707 

-.000 

-.707 

-1,000 

-,707 

,000 

,707 

•31 3 E -05 

.717 

1,053 

.772 

,039 

-.717 

-1.053 

-,772 

-.039 

.313 E -03 

,707 

1,000 

.707 

-.000 

-.707 

- I.OOO 

-.707 

.000 

•32 bt -03 

-.031 

,706 

1.029 

,749 

• 031 

-.706 

-1.029 

-,749 

•325 E -03 

• ooo 

. 707 

1,000 

,707 

-,000 

-• 707 

-1,000 

-,707 

, 336 E -03 

-.743 

-.038 

,689 

1.013 

• 743 

.038 

-,609 

-1,013 

• 336 E -^0 3 

-.707 

.000 

,707 

1,000 

,707 

-.000 

-.707 

-1,000 

•3 b 0 t -03 

- l.OlO 

-.746 

-.045 

,683 

1,010 

,746 

,045 

-.663 

.350 E -05 

-1,000 

-,707 

.000 

,707 

1,000 

.707 

• ,000 

•,707 

• 363 E -03 

- , 683 

•1,008 

-.743 

-,042 

,683 

1,006 

,743 

• 042 

•363 E -03 

-.707 

•1,000 

-.707 

,000 

,707 

1 .000 

, 7 U 7 

• ,000 

•375 E -03 

• 039 

-.677 

-.996 

-.732 

-.039 

.677 

,996 

.732 

•375 E -03 

-.000 

-.707 

-1,000 

-.707 

,000 

.707 

1,000 

,707 

• 3 B 8 E -0 3 

,7 2 b 

.047 

-.660 

-.980 

-,72 b 

-.047 

,660 

,980 

.3 B 8 E -03 

,707 

-,000 

-,707 

-1,000 

-.707 

.000 

,707 

1,000 

,<»00 E -03 

,974 

. 735 

.065 

.,643 

..974 

-.735 

-.065 

,643 

.400 E -03 

1.000 

,707 

-,000 

-.707 

-1 ,000 

-.707 

,000 

• 707 

•« l 3 E -03 

.637 

,981 

,751 

• 060 

-.637 

-.981 

-.751 

-.080 

.413 E -03 

,707 

1,000 

,707 

• ,000 

-,707 

-1 .000 

-,707 

,000 

• ^2 bE -05 

-,08 a 

,642 

.991 

,760 

,084 

-.642 

-.991 

-.760 

•425 E -03 

,000 

,707 

1.000 

,707 

-.00 0 

-,707 

-1 .000 

707 

• M 38 t -0 3 

762 

-.063 

.644 

,994 

.762 

.063 

-.644 

-,994 

•436 E -03 

-.707 

,000 

.707 

1.000 

,707 

• ,000 

-.707 

•1,000 

•« bOE -03 

-.993 

- , 764 

-.008 

,640 

,993 

,764 

,008 

-•640 

•450 E -03 

-1,000 

-,707 

,000 

,707 

1 ,000 

,707 

-,000 

• ,707 

•^63 E -03 

-.637 

-.996 

-,77 b 

-.097 

.637 

,990 

,775 

.097 

•463 fc -03 

-.707 

-1,000 

-.707 

,000 

,707 

1,000 

.707 

• ,000 

• < l 7 bE -0 5 

.102 

-,643 

•1.011 

-.787 

-,! Q 2 

,643 

1.0 M 

.767 

.475 E -03 

-.000 

707 

-1.000 

-.707 

,000 

, 707 

1,000 

• 707 

, <<60 E - O 3 

.791 

,09 b 

-.657 

-1,024 

-.791 

-.095 

,657 

1,024 

•408 E - O 3 

,707 

-.000 

-.707 

-1,000 

-.707 

• 000 

,707 

1,000 

• bOOt -03 

1.020 

,767 

,085 

.,667 

. 1.020 

-.787 

-,005 

.66/ 

. bOOE-Oi 

i.ooo 

.707 

-.000 

-.707 

•1,000 

-.707 

• 000 

,707 



TIME 

GREENS 

FUNCTION HR 



TIME 

GREENS 

Function hh 



0.000 

25,065 

-.000 

-,000 

-,ooo 

0,000 

8. on 

-,003 

-.056 

,075 

,00a 

-.005 

.157 

-1,681 

,070 

,260 

.317 

.320 

6,160 

-.001 

-.037 

,069 

,009 

,000 

.310 

-2.305 

.165 

.370 

.596 

.000 

8,325 

-.070 

-.030 

.055 

,010 

.003 

.071 

-2,630 

.332 

,389 

.395 

,397 

6.062 

-.060 

-.010 

,038 

,006 

.002 

.6ae 

-2,608 

,SOH 

,301 

,518 

,310 

8,639 

-.055 

,006 

,021 

,002 

-.003 

,7dS 

•2,000 

.078 

,200 

,190 

,185 

6,796 

-.052 

,027 

,009 

-.005 

-,009 

.9oa 

-1,999 

,800 

,m 

,060 

,050 

8,950 

-.057 

,005 

,001 

-.011 

-,oia 

1,100 

-1.525 

,055 

-.016 

-.009 

-.055 

9,m 

- , 067 

,056 

-.001 

•.013 

-.016 

1,2S7 

-1,001 

.015 

-.113 

-,lOb 

-.107 

9,260 

-.079 

,059 

,000 

-.01? 

-,015 

1,010 

-.731 

,000 

-.159 

-.111 

-,100 

9.025 

-.007 

,055 

.003 

-,00d 

-•on 

1.S71 

-.505 

,097 

-.106 

-.069 

-.050 

9,582 

-.005 

• 006 

,005 

-.005 

-.008 

1.726 

-.395 

,200 

-.077 

-,007 

.003 

9.739 

-.071 

,035 

,006 

• ,003 

-,005 

1.805 

-.367 

,076 

.050 

,051 

,050 

9,896 

-,0OO 

,020 

,oou 

-,003 

-.005 

2,002 

-.379 

-.083 

,153 

,082 

,076 

10,053 

-.009 

,015 

- , 000 

• ,006 

-.008 

2.199 

-.392 

-.180 

,267 

,077 

.060 

10,210 

.025 

,009 

-.005 

-,010 

-,012 

2,356 

-.300 

-.215 

.351 

,000 

,020 

10.367 

.053 

.003 

-.008 

-,oia 

-.015 

2,513 

-.300 

-.197 

,390 

-.011 

-.016 

10,520 

.067 

-.000 

-.010 

-.015 

•,016 

2,670 

-.293 

-.155 

,382 

-.056 

-,007 

10,601 

,06b 

-.003 

-.010 

•,01a 

-.015 

2.827 

-.236 

1 1 0 

.332 

-.076 

-.050 

10,630 

.052 

-.003 

-,007 

-,010 

-.012 

2,985 

-.191 

-.078 

.255 

-.060 

-.030 

10,996 

,029 

,000 

-,004 

-.007 

-,0U8 

3.102 

-.167 

• , 067 

.167 

-.018 

,005 

11.153 

.005 

,007 

-.002 

-.005 

- ,006 

3,299 

-.162 

-.075 

.065 

,052 

,039 

11,310 

-.015 

,017 

-.001 

-,005 

-.006 

3.056 

-.170 

-.091 

.020 

.126 

,059 

11,067 

-.011 

,001 

-,003 

-,020 

• ,026 

3.613 

-.179 

-.106 

-.020 

.189 

,05 3 

11,620 

-.021 

,000 

-,007 

-.019 

-,02a 

3,770 

-.101 

-,lll 

-.037 

,226 

,021 

11.761 

-,02b 

.052 

-.Oil 

-.018 

-,020 

3.927 

-.172 

-.103 

• . 036 

.237 

-.023 

11,950 

-,02b 

,053 

-.015 

-,016 

-,017 

0.060 

-.153 

-.087 

-.025 

,219 

-.061 

12.095 

-, 020 

.050 

-.015 

-.015 

-.016 

0,201 

-.132 

- ,067 

-.013 

,103 

-.075 

12,252 

-.021 

,005 

-.012 

-,oia 

- , 0 16 

0,598 

-.110 

-,051 

- , 006 

.159 

-.051 

12,009 

-.020 

,G36 

-,005 

-,oia 

-.018 

0,555 

-.100 

-.000 

-.006 

,097 

.011 

12,566 

-.021 

.031 

.005 

-.015 

-.021 

0,712 

-,ioo 

-.005 

-.010 

,061 

.100 

12,723 

-.020 

,020 

,016 

-,017 

-,0?a 

0,869 

-.no 

-,052 

-.026 

,050 

,197 

12,881 

-.029 

,019 

,027 

-.019 

-.025 

5,027 

-.117 

-.060 

-.057 

.015 

,280 

13.038 

-.030 

,015 

,03a 

-,020 

- , 0?5 

5,180 

-.120 

-,060 

-,003 

,00? 

.529 

13.195 

-.037 

.013 

.038 

-.021 

-.023 

5,301 

-, 1 16 

-.060 

-,005 

-.007 

,336 

13.352 

-.023 

,013 

,021 

-.020 

-.ooa 

5,090 

-.106 

-.051 

-.036 

-,Oll 

,503 

15,509 

-.027 

.015 

,020 

-.016 

-.003 

5,655 

-.093 

-.000 

-.027 

-.0 08 

.20? 

13,666 

-.030 

.013 

,017 

-.013 

-,005 

5.812 

-.083 

-.030 

-.017 

,000 

,170 

13.623 

-.032 

,015 

,013 

-,007 

-.000 

5,969 

-.0 76 

-.025 

-,012 

,0?0 

.100 

13,960 

-.053 

.013 

,006 

,001 

-.012 

6.126 

-.079 

- ,026 

-.01? 

,050 

,056 

10,137 

-.032 

.012 

.003 

,008 

-.016 

6,203 

-.005 

-.032 

-.018 

.077 

.031 

10.290 

-.031 

• Oil 

-.001 

.015 

-.019 

6,000 

-.091 

-.036 

-.020 

.099 

,025 

10,051 

-.029 

,009 

-.005 

,020 

-.020 

6,597 

-.090 

-.00? 

-.037 

• 111 

.030 

10,606 

-.026 

,008 

-.000 

,02? 

-.018 

6,750 

-.092 

-,0Ol 

-,0U2 

,110 

,037 

10.765 

-.027 

.007 

-.009 

.022 

-.012 

6,912 

-.000 

-,0 3« 

-.039 

,095 

.059 

10.923 

-.027 

,006 

-,oin 

.019 

-.003 

7,069 

-,075 

-.025 

-.027 

,07? 

,030 

15,060 

-.023 

.006 

-.009 

,Oia 

,007 

7,226 

-.067 

-.016 

-.008 

,Oub 

,022 

15,237 

-.050 

• OOt 

-.008 

,007 

,020 

7,303 

-.060 

-.011 

,017 

,02? 

,000 

15,390 

-.032 

,007 

-.007 

.001 

,031 

7,500 

-.066 

-, 0 1 3 

.001 

,006 

-.000 

15.551 

-.030 

,007 

-.006 

-.005 

.0U2 

7,697 

-.072 

-.020 

.061 

-,001 

-.010 

15.706 

-.035 

.007 

006 

-,009 

,oa9 

7,650 

-.079 

-.029 

,073 

• , 000 

-.009 








KjJ 

"vj 


\JV> 

OO 


TIME 

GHEENS 

FUNCTION 

TIME 

GREENS 

FUNCTUiN JR 

0.000 

•3,698 

-1,520 

-.629 

8, on 

,016 

-.129 

-.033 

.157 

-3,114 

•1.248 

-.514 

6,166 

,024 

-.116 

-.033 

. il 4 

-2.209 

-.851 

-.346 

8,32 b 

,023 

-,093 

-,029 

.471 

-1.138 

-.420 

-.171 

0.462 

,011 

-.066 

-.025 

.626 

-,07 b 

-.054 

-.024 

6.639 

•,0 J 5 

-,044 

-.017 

,765 

.620 

,165 

,072 

8,796 

-.045 

-.032 

-.013 

,9 a 2 

1,446 

,264 

,098 

6,954 

-,077 

-,029 

-,012 

1,100 

1,764 

,196 

,073 

9. in 

-.101 

-.032 

-.014 

1,257 

1.788 

,050 

,020 

9,268 

-.113 

-.037 

-,016 

l , ai 4 

1,590 

-.106 

-,05? 

9,425 

-.109 

-.039 

-,016 

1,571 

1,267 

-.208 

-,059 

9.582 

-.093 

-.037 

-• 0 16 

1,728 

.912 

-.203 

-.055 

9,739 

-,069 

-.030 

-.013 

1.66 S 

,601 

-.090 

-,024 

9,696 

-,046 

-.022 

-.010 

2,042 

.375 

,105 

,017 

10,053 

-.029 

-,016 

-.008 

2,199 

.242 

.329 

,046 

10.210 

-,023 

-.014 

-.007 

2,356 

.184 

,525 

,050 

10,367 

-.027 

-.016 

-.007 

2,513 

.172 

,650 

,024 

10,524 

-.037 

-,020 

-.009 

2,670 

.1 7« 

,601 

-,020 

10,661 

-.046 

-,023 

-,OlO 

2,627 

• 170 

,625 

-.062 

10,636 

-.049 

-.024 

-,on 

2,965 

,152 

,510 

-.080 

10,996 

-.041 

-.020 

-.009 

^,142 

» I 23 

.373 

• • 056 

11,153 

-.024 

-,015 

-.007 

3,299 

.091 

,249 

,012 

11,310 

,000 

-,009 

-.005 

3,456 

,065 

,163 

,114 

11,467 

,019 

-.015 

-.010 

3,613 

,051 

• 121 

,226 

11,624 

,057 

-.020 

-.013 

3,770 

,048 

.113 

,325 

11,761 

.050 

-.025 

-,013 

3,027 

,052 

,122 

,383 

11,938 

,057 

-.023 

-.oil 

4,064 

,055 

.131 

.390 

12,095 

,056 

-.020 

• ,006 

4,241 

.054 

,126 

.340 

12,252 

,048 

-.012 

-.001 

4,396 

,046 

• Ul 

,273 

12,409 

.037 

-.002 

,002 

4,555 

,037 

,085 

,168 

12,566 

,025 

• on 

,003 

4,712 

,027 

,059 

• in 

12,723 

,015 

,024 

,001 

4,669 

.020 

,043 

,075 

12.661 

,006 

,036 

-.004 

5,027 

,019 

,039 

,066 

13,038 

,006 

,045 

-.010 

5,184 

• 021 

,046 

,082 

13,195 

.007 

,050 

-.015 

5,341 

.025 

,055 

,108 

13.352 

,004 

,051 

-,010 

5.496 

.027 

,061 

,125 

13,509 

.004 

.046 

-.003 

5,655 

.026 

,056 

,120 

13.666 

,005 

.042 

,006 

5,612 

,021 

,045 

,069 

13.623 

,008 

,034 

.016 

5,969 

• 014 

,029 

.036 

13,960 

• Oil 

,026 

.025 

6,126 

,009 

,016 

-.027 

14,137 

,013 

,019 

,033 

6,263 

• 007 

,011 

-,086 

14,294 

,014 

,014 

.039 

6,440 

,009 

,017 

-.127 

14.451 

.013 

• on 

.042 

6,597 

,013 

,031 

-.145 

14,606 

.on 

,011 

,042 

6, 754 

,016 

,045 

-.139 

14,765 

,006 

,012 

.041 

6,912 

,019 

,051 

-.116 

14,923 

,005 

,014 

,039 

7,069 

,015 

,042 

-.066 

15.060 

.003 

.016 

,035 

7.226 

,009 

,016 

-.057 

15,237 

,001 

,017 

,032 

7.363 

,003 

-,025 

-.037 

15.594 

,001 

,016 

.028 

7.540 

-.001 

-,066 

-.028 

15,551 

.001 

.015 

,025 

7.697 

,001 

-,103 

-.027 

15.706 

,003 

,012 

.022 

7,654 

,008 

-,125 

-,030 







TIME 

GREENS 

EU-'iClIUN KR 



TIME 

GREENS 

FUNCTION KW 



0.000 

25,465 

•,000 

-.000 

-.000 

0,000 

6,011 

-.009 

-.004 

,008 

,000 

-.001 

• 157 

•,167 

,009 

,052 

,055 

,056 

6.166 

-.009 

-,004 

,003 

,001 

,000 

• JH 

•, 256 

,020 

,041 

,044 

.045 

6,325 

-.008 

-,00 5 

,006 

,001 

,000 

• <<7l 

-.292 

.057 

,043 

.044 

,044 

6,482 

-.007 

-.002 

.004 

,001 

,000 

, b 2 S 

-.294 

,056 

,036 

.055 

.035 

6,659 

-,00b 

.001 

,002 

,000 

-.000 

, 765 

•,266 

,075 

,027 

,022 

,021 

6,796 

-.006 

,003 

.001 

-,C01 

-.001 

.942 

•.222 

,069 

,012 

,007 

,006 

6,954 

-.006 

,005 

,000 

-.001 

- ,002 

1.100 

-.169 

,095 

-.002 

-.005 

-.006 

9,111 

-.007 

,006 

-.000 

-,00l 

-.002 

1.257 

•.120 

,091 

-,013 

-.012 

-.012 

9.266 

-.009 

,007 

.000 

-.001 

-.002 

l.ai4 

-,08l 

,076 

-.018 

-.012 

-.012 

9,425 

-.010 

,006 

.000 

-.001 

-,001 

1.571 

•,056 

,055 

-,016 

-,006 

-.006 

9,562 

-.009 

,005 

.001 

• , 0 0 1 

-.001 

r.726 

-.044 

.051 

-.009 

-,0U1 

,000 

9,759 

-.006 

,004 

.001 

-,00O 

-.001 

1,665 

••041 

,006 

,003 

,006 

,006 

9,696 

-.005 

,005 

,000 

-.00 0 

-.001 

2,042 

-.042 

-.009 

.017 

,009 

,008 

10.053 

-.001 

,002 

-,000 

-.001 

-,001 

2,199 

•,044 

-.020 

,030 

,009 

,007 

10,210 

.005 

,001 

-.001 

-,001 

-,001 

2,i56 

-.043 

-.024 

.059 

,004 

.003 

10.567 

.006 

,000 

-.001 

-.002 

-.002 

2,515 

-.059 

-.022 

,043 

-,001 

-.002 

10,524 

.007 

-.000 

-,001 

-.002 

-.002 

2,670 

-,053 

-,017 

',042 

-,006 

-.005 

10.661 

.007 

-.000 

-.001 

-.002 

-.002 

2,627 

-.026 

-.012 

,057 

-,009 

- , 006 

10,838 

,006 

-.000 

-,001 

-,001 

-,001 

2,965 

-.021 

-.009 

.028 

-.007 

-.004 

10,996 

,003 

.000 

-.000 

-.001 

-,ooi 

5,142 

-.019 

-.007 

.019 

-,002 

,000 

11.153 

.001 

.001 

-,000 

-.001 

-.001 

5,299 

-.018 

000 

,009 

,00b 

,004 

11.310 

-.002 

.002 

-.000 

-.001 

-.001 

5,456 

-.019 

-.010 

,002 

,014 

,007 

11,467 

-.001 

,005 

-,000 

-.002 

-.005 

5.615 

-.020 

-.012 

-.002 

,021 

,006 

11,624 

- , 002 

,005 

-,001 

-.002 

-,005 

3,770 

-.020 

-.012 

-.004 

,025 

.002 

11.781 

-.005 

,006 

-,001 

-.00? 

-.002 

5,927 

•,019 

".oil 

-,C04 

,C2b 

-,005 

11.936 

-.005 

,006 

-.002 

-.002 

-,002 

4,0fi4 

-.017 

-.010 

-.005 

,024 

-.007 

12,095 

-.005 

,006 

-.002 

-,002 

-,002 

4,241 

-.015 

-.007 

-.001 

,020 

-.008 

12,252 

-.0 02 

.005 

-.001 

-.002 

-.00? 

4,596 

-.015 

•,00b 

-.001 

.015 

-.006 

12,409 

-.002 

,004 

-,001 

-.002 

-,0l)2 

4,555 

-.012 

• ,005 

-.001 

,011 

,001 

12,566 

- , 002 

.005 

.001 

-.002 

-.002 

4,712 

-.012 

-.005 

-.002 

,007 

,011 

12,723 

-.003 

• 00 3 

,002 

-,0U2 

-.005 

4,669 

-.012 

• , 006 

-.00 5 

,004 

,022 

12,861 

-.003 

,002. 

.005 

-,002 

-,003 

5,027 

-,015 

-.007 

-.004 

,00? 

.051 

15,038 

-.004 

,002 

,004 

-.002 

-.003 

5,164 

-.015 

-.007 

- , 005 

,000 

,037 

13,195 

-.004 

.001 

,004 

-.002 

-.003 

5,341 

-.013 

-.007 

-,005 

-.001 

,037 

13,352 

- , 005 

.001 

.002 

-,002 

-.000 

5,496 

-.012 

- , OOh 

-,00a 

-.001 

,034 

13.509 

-.003 

,001 

.002 

-.002 

-.000 

5,655 

-.010 

-,OOU 

-.005 

-.001 

,027 

13,666 

- , 005 

.001 

,002 

-.001 

-.001 

5,812 

-.009 

-.005 

-.002 

,000 

,019 

13,823 

-.004 

.001 

,001 

-.001 

-.001 

5,969 

-.009 

-.005 

-.001 

,003 

.012 

13,960 

-.004 

.001 

,001 

,000 

-.001 

6,126 

-.009 

•,005 

-.001 

,006 

.006 

14,137 

-.004 

,001 

,000 

,001 

-.002 

6,265 

-,0U9 

-.004 

-.00? 

,009 

,003 

14.294 

- , 005 

,001 

-.000 

,002 

-,002 

6,440 

-.010 

-.004 

-.005 

,011 

,003 

14,451 

-.005 

.001 

--.001 

,002 

-.002 

6,597 

-.010 

-,005 

- ,004 

.012 

,003 

14,603 

-.005 

.001 

-.001 

,002 

-.002 

6.754 

U 1 0 

-.005 

-.0 05 

,012 

,004 

14.765 

-.005 

,001 

-,00l 

,002 

-.001 

6.912 

-.009 

-.004 

- , f ’ 0 4 

.011 

,004 

14.925 

-.003 

,001 

-,001 

,002 

-,000 

7,069 

-.006 

-.005 

-.003 

,008 

,0 04 

15,060 

- , 005 

.001 

-.001 

.002 

,001 

7,226 

-.007 

-.002 

-.001 

,0 05 

,002 

15.237 

-.005 

.001 

-.001 

,001 

,002 

7.365 

-.00? 

• , 00 1 

,002 

,002 

,001 

15. 394 

-.004 

,001 

-.001 

,000 

,005 

7,540 

-.007 

-,coi 

,0 05 

.001 

-,000 

15,551 

-.004 

,001 

-.001 

-,001 

,005 

7,697 

-.006 

-.002 

,00 7 

- , 0 0 0 

-,001 

15,708 

-.004 

,001 

-.001 

-.001 

• oos 

7,654 

-.009 

-,C05 

,008 

-.000 

-.001 








VjO 

\ S > 



TIME 

PRESSURE P(TIME) 

FOR DISCRETE VALUES 

OF Y 




0, 


.526EtOS 

,372E+06 

-.760E-09 

-.3726+06 

-.526E+06 

-.372E+06 

-.822E-06 

.3726+06 


0« 

.427E+06 

.570E+06 

.379E+06 

-,33^E+05 

-.427E+06 

-.570E+06 

-.5796+06 

.3396+05 

• asoc 

0« 

.2211+06 

,603E+0b 

. 63 1 E + Ob 

,2906+06 

-.221E+06 

-.603E+06 

-.651 E+06 

-.2906+06 

• 573t 

OJ 

2out + 06 

•254E+06 

.624E+06 

.6286+06 

.2646+06 

-.2546+06 

-.6246+06 

-.6286+06 

.soot 

0 4 

-.6R2E+U6 

-.521E+06 

,236L*0b 

,658t+06 

.6926+06 

.321E+C6 

-,258£+P6 

-.6586+06 


0 4 

“■62it+ 06 

-.843E+06 

“,368E+06 

.3226+06 

,8256+06 

.6436+06 

.568E+06 

-.3226+06 

.7S0E 

OJ 

-,550L*nb 

-.102E+07 

-,909t*06 

-,268t+06 

,5506+06 

,1026+07 

.909E+06 

, 2b8E + 06 

.67SE 

04 

.296Etu5 

-.752E+06 

-.1 09E + 07 

•*,794Et06 

-.296E+05 

.752E+06 

.1096+07 

, 7946+06 

.loot 

05 

,S76t+06 

1 76L + 06 

-.825E+06 

-,99lt+06 

-,5766+06 

, 1 76E + 06 

,8256+06 

.99U + 06 

.I13t 

05 

.0O2£*Ob 

,597E+06 

-.2U1E+06 

-.7376+06 

-.8026+06 

-.597E+06 

,2416+06 

.7376+06 

. ust 

05 

.609E+06 

,679E+Ob 

.351 t+06 

-.1 82t + 0b 

6096 + 06 

-.679E+06 

-.351E+06 

,1826+06 

,lidt 

05 

, 109£*06 

,557Et06 

,651 1 + 06 

.363E+06 

-.1096+06 

-.557E+06 

-.651E+06 

-.3936+06 

.ISOt 

05 

-.3'73£tC6 

.971E+0S 

.53CE+06 

,6556+06 

.3956+06 

-.9716+05 

-.5506+06 

-.6536+06 

« 1 bit 

05 

-.620E+C6 

•,5b5t»fi6 

. 1 05E + 06 

,5136+06 

,6206+06 

.565E+06 

-.1056+06 

-.5136+06 


05 

•, >*4t;E + 0b 

-,559£f 06 

-.3U4E+06 

,7216+05 

,4466+06 

.559E+06 

, 34UE + 06 

-.7216+05 

. 10bt 

05 

,7l9E*0a 

-,375E+06 

-.534E+06 

-,3836+06 

-.7196+04 

,5756+06 

,5346+06 

.3956+06 

• ^ 0 0 L 

03 

.4bOE*06 

•797E+05 

-.347E+06 

-.570F+06 

-.460E+06 

-.797E+05 

, 5476+06 

,5706+06 

,^lit 

05 

, b26£ ♦ 06 

.515E+06 

.lOOE+06 

-.3736+06 

-,628E-* 06 

-.5156+06 

-.1006+06 

,3736*06 

• 2i5£ 

05 

.40u£+06 

, 6 5 9 E ♦ 0 6 

.527E+06 

,8b8E+05 

-.404E+06 

-.6596+06 

-, 5276 + 06 

-.9686+05 

.^iat 

05 

-,{Ju7£tC5 

,4C8t+06 

.661E+06 

,5276+06 

,6476+05 

-.408E+06 

-.6616+06 

-.5276+06 

. 2Sot 

05 

•, SS 1 £ + 06 

-,lP5t+06 

,402t+06 

.6736+06 

.551E+06 

. 1 05E + 06 

-,u02t*06 

-,67it+06 

,^bit 

0 5 

•, 7C6E + C6 

-.581 E+06 

1 15E+06 

.4186+06 

,7066+06 

,5816+06 

, 1156 + 06 

-.uibt+Ob 

,27!:t 

05 

•.457E+C6 

-.742E+06 

-.595E406 

-.9836+05 

,4S7fc+06 

.7426+06 

.5956+06 

.9636+05 

.2r6t 

05 

,S7o£+05 

• . 49 5E +06 

-,755t+06 

-.5756+06 

-.5766+05 

.495E+06 

.7556+03 

.5756+06 

. 300E 

05 

,542c ♦ 06 

.268E+05 

-.504E+U6 

•,7U0F +06 

-.5426+06 

-.268E+05 

,5046+06 

,7406+06 

,513l 

05 

.718E+06 

.521E+06 

.188E+05 

-.4946+06 

-.7186+06 

-.5216+06 

-.1886+05 

. u9iit*0b 

.32St 

05 

.aoSE+Ob 

.7O0E+O6 

.517E+06 

,2286^^05 

-.4856+06 

-.708E+06 

-.5176+06 

-.2296+05 

.35bE 

03 

-.205E+05 

,486fc+06 

.708E+06 

,5156+06 

.2056+05 

-.466E+06 

-.708E+06 

-.5156+06 

. isot 

05 

-.504E+06 

•.990E+04 

.4906+06 

.702E+06 

.5G4E+06 

,990E+04 

-.4906+06 

-.7026+06 

,36it 

05 

b85E + 06 

-.4B7E+06 

•,449t+04 

.481E+06 

,6656+06 

.4976+06 

,4496 +U4 

-.4816+06 

.37St 

05 

•,462E+0o 

667E+06 

•,461 t+06 

-,1346+05 

,4626+06 

.667E+06 

,4816+06 

.1346*05 

.3aat 

05 

,2<?b£*C5 

•,447e+U6 

•.662t+06 

- ,4896 + 06 

-.2966+05 

,4476+06 

,6626+06 

.4896+06 

.uoot 

05 

,a99£tOb 

.398E+05 

-,4u3c+06 

-,6666+06 

-.4996+06 

-.396E+05 

,4436+06 

.6666+06 

.Ulit 

05 

, b 7 l) t ♦ 0 0 

.S05E+06 

.416E+05 

-.44UE+06 

6706 + 06 

-.5056+06 

-.4156+05 

, 4446 + 06 

.4^5£ 

03 

, 4 4 1 E ♦ 0 0 

, 6b7£ ♦ 06 

.502E+06 

,435E+05 

-,U4lE+0b 

-.6676+06 

-.5026+06 

-.4356+05 

.43bE 

05 

-,526Et05 

,453E»06 

.66UE+06 

,S07E+0b 

.525E+05 

4i3E + 0b 

6646 + 06 

-.5076+06 

.4O0E 

05 

-,519£+06 

-,b57tt05 

.U29£i06 

,6706+06 

.5196+06 

.637E+05 

-.4296+06 

-.6706+06 

• «t>3E 

05 

•, 68 5E + 06 

•,53 lE *06 

-.6766+05 

.4356+06 

.665E+06 

,531£+06 

,6786*05 

-.4356+06 

.4?S£ 

05 

•,445E+06 

-.692E+06 

-.534E+06 

-.6326+05 

,4456+06 

.6926+06 

,5346+06 

. 6326+05 

,«ddE 

05 

.572E+05 

•.asit+06 

-.695E+06 

-.532E+06 

-.5726+05 

.451E+06 

• 695E+0b 

.5326+06 

• SOOt 

03 

.551E+06 

.5626+05 

-,o5lt+06 

-.69UE+06 

-.5516+06 

-.5626+05 

. 4516 + 06 

.6946+06 



time 

UPSTREAM COMPONENT OP PRESSURE 





0. 


,500Et06 

,354E*0h 

-.858E-09 

-■354E+06 

-.500E+06 

-,35ilE + 06 

•,e09E-08 

.3SaE+06 

.USE 

04 

,327E^0b 

,475E*06 

.345E+06 

.l23Et05 

-.327E+06 

-.475E+06 

-.345E+06 

-,U3t + 05 

.250E 

00 

•.IIRE^OS 

.327E*06 

.47tiE + 06 

.543E+06 

,119E+05 

-.327E+06 

•,47UE+06 

-.343E+06 

.375E 

oo 

-,349E+0b 

-.167E+05 

.326E+06 

.a77fct06 

.34QE+06 

.lb7t+05 

-.3261+06 

-,477t+06 

.500E 

04 

-.uSOEtOb 

-.352E+06 

-.160E+05 

,327t406 

,«80E+06 

.3S2E+06 

.180E+05 

-,327E+06 

,02S£ 

04 

-.5531*06 

-,46fct+06 

-.3SUE+06 

-.150E+05 

.335E+06 

,4B6E-+06 

,3S4E+06 

.150E+05 


0** 

.iUE*0S 

-.53bt*06 

-.4e7t*06 

-.352E+06 

-.112t+C5 

.336E+06 

,487E+0b 

.352E+06 


Oa 

,347Et0b 

,650i*0-4 

•.35eE*06 

-.4flSE+06 

-.3U7E+06 

-.6S0E+04 

.338E+06 

,U8St+06 

.loot 

03 

,4S2E*0b 

,3«4E+06 

.547E+04 

-,357E+06 

-,4B2E+06 

-.344E+06 

-.5U7E+04 

.337E+06 

,113£ 

03 

.33«E*0b 

.479E+06 

,343E*0b 

.691E+04 

-.334E+06 

-.479E+06 

-.343E+06 

^,69lE+04 

.USE 

05 

-.80SE+04 

.333E+06 

.476E+06 

.344E+06 

,8081+04 

-.333E+06 

-,476E+0b 

-,34ttt+06 

,136t 

03 

-.5451*06 

-.916E+04 

.532E+06 

.U79E+06 

,345t+06 

, 9181 + 04 

-,332t+06 

-.479fc+0b 

.iSuE 

03 

-.a7REt0b 

-,34SE*06 

-.909E+04 

,352E+06 

,u79e +06 

.345E+06 

.9O9E+04 

-,332E+06 

.UiE 

03 

-,332E*06 

-.478t*06 

-.34SE+06 

-.922t+ua 

,3321+06 

,U78£+06 

.345E+06 

.922E+04 

. wst 

03 

.loot *0S 

-,331t*06 

-.47PE+06 

-.34SE+06 

-.lOOE+05 

.33U + 06 

,U78fc+06 

,3ubE+06 

.USE 

0 5 

,346E:*06 

•107E+05 

-, 5311 + 06 

-.479E+06 

-,34bt+06 

-.1071+05 

,331E+06 

, 4 7bt + 06 

.200 1 

03 

, 479E + 06 

,3471+06 

.109E+0S 

-.33U + 06 

-,u79t+06 

-.3U7E+06 

-.109E+05 

,3ilt+06 

.2lit 

03 

.5521*00 

.aeot *06 

,5«7E*06 

, lOflt +05 

-.332E+06 

-.UBOE+06 

-.347E+06 

-.106E+05 

.22St 

05 

-,107E*0S 

,332t*06 

,460E*06 

,347t+06 

.107t+0b 

-,332E+06 

-,4e0E+06 

-.347E+06 

.23SE 

03 

-.3«7E*06 

-.1 10£*0S 

.352E+06 

, u80E +06 

.347E+06 

.llOE+OS 

-,332t+0b 

-,u60E+06 

.2S0E 

03 

-,4eyE*0o 

•,54riF+0b 

-.lUt + OS 

.332E+06 

,U0OE+O6 

,34BE+06 

.lUE + 05 

-.352E+06 

,ibSt 

03 

-,352E*0b 

•,UE0t*06 

-,34et*06 

-,n2£ + 05 

.332E+06 

,u60E+06 

.348E+06 

.112E+0S 

.27SE 

05 

.lUE + 05 

-.3321+06 

-.480E+06 

-, 5U8E + 06 

-.112E+05 

.332E+06 

.490E+06 

,3yeE+06 

.2SSE 

05 

,3461*06 

,112l*05 

-,332E‘+06 

-.U60E+06 

-,34Pt+06 

-.112E+05 

,3521+06 

,490t+06 

. 30oE 

05 

.uiJOEfOo 

, it) t E ♦ 06 

• 1 1 2t + CS 

-.332E+06 

-.480E+06 

-,348E+06 

-,M2E + 0S 

.332E+06 

,51 it 

03 

.352E*06 

,460E*Cb 

,346t*06 

.M2E + 05 

-,3*32t+06 

-.UBOE+Ob 

-.348E+06 

1 1 2E + 05 

.32SE 

C3 

1 UE*05 

,3321+06 

. 4d0£ + 06 

, 5udt + 06 

.112E+0S 

-.332E+06 

-,480E+06 

346E+06 

.33St 

05 

346c *06 

-,112t+0S 

. 552E+06 

.460E+06 

,346f +06 

.lUE + 05 

-.352E+06 

-,u60£+06 

.3S0E 

05 

- , a 6 0 1 ♦ 0 6 

-. 34CF + C6 

-.1 12E + 0S 

.352E+06 

.4B01+06 

.346E+06 

• lUE + OS 

-,532E+06 

• ioit 

05 

-,352E*06 

-,u60E*Ob 

34Pt + 0b 

-.lUfc+OS 

.552t+06 

,40Ot+Ob 

,3u6E+06 

,112E+0S 

• 37St 

0 5 

.112t*05 

-,532E+C6 

-.4P0E+06 

-,346E+06 

-,U2t + 0S 

,332t+06 

,4b0E+06 

.346E+06 

.30{?t 

03 

, 34 9E ♦ Ob 

.112E+C5 

-.3321+06 

-,4BnE+06 

-.340E+O6 

-.112E+05 

,332t+06 

,4bUE+06 

,**OOE 

03 

.46. It *06 

, i u t E * 0 6 

.1 Ufc + OS 

-.332E+06 

-.4901+06 

-,34dE+06 

-.lUt + 05 

, 352E + nt, 

,yl <t 

C5 

,332£*0b 

,4.5 0t*0b 

,54'’E + 0b 

.lUt + nS 

-.3521+06 

4b0 1 + 06 

-.348E+06 

-.112E+05 


03 

-,112£*05 

.332E+C6 

,4601+06 

,3481+06 

.1 12E+0S 

-.332E+06 

-,4fl0t+06 

3U6E + 06 

,u5St 

C5 

5USE • 06 

•,112h+05 

.332E+06 

,4P0t+06 

, 5UBL + 06 

. 1 12E + 05 

-.352E+0b 

-.48CE+06 

• *iS0t 

05 

-. 460£ *06 

-,3a3t*06 

H2t + 05 

.5521+06 

.U60E+06 

, 348E+06 

.lUc + OS 

-,i52E*06 

.Obit 

03 

-.332E*0b 

-.JeCE+Cb 

3ubt + C6 

-.112E+05 

,332£+06 

.4B0E+06 

, 346E+06 

,n2fc + 05 

.M /bL 

03 

, 1 Ut*0b 

-,i52£+06 

-.4601+06 

-,34UE+0b 

-.M2E + 05 

.332E+06 

.480E+06 

, 346E+06 

.««bt 

03 

, 34SE + 06 

.ll?E*0b 

-,332E+06 

-.4POE+06 

-.348E+06 

-.lUt + OS 

,332£+0b 

.ueuE+Ob 

• soot 

03 

.u60E*0b 

.3461+06 

,1 12E + 0S 

-.532E+06 

-.4601+06 

-.346E+C6 

-,M2t*05 

.352E+06 


PARAMETERS 

USED 

OTlMt 

E 

.250E-05 

NY 

E 

6 

NT 

= 

201 

MAXCON 

S 

lOl 

CMAKK 

S 

.eooE+00 

RHOO 

s 

, looe^oi 

C 

B 

. looE+oa 

Y 

3 

. lOOE^OO 

OMEGAO 

X 

,S00E400 


TEST RISTAR approaches ITS ASYMPTOTIC VALUE 


TIME 

STEP 

1 

RISTAR 

3 0.000 

,1413 

ASYMPTOTE 

3 

2,500 

1,658 

time 

STEP 

6 

RISTAR 

* *.166 

.827 

ASYMPTOTE 

B 

1.675 

2,489 

TIME 

STEP 

11 

RISTAR 

3 -.575 

1,355 

ASYMPTOTE 

= 

,595 

2,940 

TIME 

STEP 

16 

RISTAR 

= -I.l«3 

1 ,608 

ASYMPTOTE 

= 

-.575 

2,944 

TIME 

STEP 

21 

RISTAR 

* •1,750 

1,526 

ASYMPTOTE 

r 

-1,658 

2,500 

time 

STEP 

2b 

RISTAR 

a -2.262 

1.102 

asymptote 

3 

•2,489 

1.675 

TIME 

STEP 

3l 

RISTAR 

3 -2.5541 

.397 

ASYMPTOTE 

3 

-2,940 

.595 

TIME 

STEP 

36 

RISTAR 

* -2.S37 

-.475 

asymptote 

3 

-2,944 

-,575 

TIME 

STEP 

*41 

RISTAR 

» -2.171 

-1.369 

asymptote 

3 

•2.500 

-1.656 

TIME 

STEP 

4lb 

RISTAR 

3 -1.4J79 

-2.12*4 

ASYMPTOTE 

3 

-1.675 

-2.489 

time 

STEP 

51 

RISTAR 

- -.593 

-2,602 

asymptote 

= 

-,595 

-2.940 

TIME 

STEP 

56 

RISTAR 

3 ,508 

-2.703 

asymptote 

C 

.575 

-2.944 

TIME 

STEP 

61 

RISTAR 

3 1.516 

-2,390 

asymptote 

3 

1.656 

-2,500 

time 

STEP 

66 

RISTAR 

a 2.529 

-1.695 

asymptote 

3 

2,489 

-1,675 

time 

STEP 

71 

Hi'.TAR 

3 2.807 

-.714 

asymptote 

3 

2,940 

-,595 

TIME 

STEP 

7b 

RISTAR 

3 2,669 

,406 

asymptote 

B 

2.944 

,575 

TIME 

STEP 

ai 

RISTAR 

3 2,a93 

1,4492 

ASYMPTOTE 

B 

2,500 

1,658 

TIME 

STEP 

66 

RISTAR 

3 1,727 

2,369 

asymptote 

3 

1.675 

2,489 

TIME 

STEP 

91 

RISTAR 

= ,662 

2,894 

asymptote 

3 

.595 

2.940 

TIME 

STEP 

96 

RISTAR 

3 -.U84» 

2.977 

asymptote 

3 

-.575 

2,944 

TIME 

STEP 

101 

RISTAR 

3 -1,592 

2.597 

asymptote 

3 

-1,658 

2.500 

TIME 

STEP 

106 

RISTAR 

a -2,q66 

1,807 

asymptote 

= 

-2,489 

1.675 

time 

STEP 

111 

RISTAR 

a -2,972 

.725 

ASYMPTOTE 

3 

-2,940 

.595 

time 

STEP 

116 

RISTAR 

* -3.021 

-.463 

ASYMPTOTE 

3 

-2,944 

.,575 

TIME 

STEP 

I2l 

RISTAR 

a -2,603 

-1.626 

ASYMPTOTE 

3 

•2,500 

-1.658 

TIME 

STtP 

126 

RIST AR 

a -1,778 

-2,527 

asymptote 

3 

-1,675 

-2.489 

TIME 

STEP 

131 

RISTAR 

3 -.672 

-3,040 

ASYMPTOTE 

3 

-.595 

-2,940 

TIME 

STEP 

136 

RISTAR 

a ,5^6 

-3,083 

asymptote 

3 

,575 

-2,944 

TIME 

STEP 

1*41 

RISTAR 

3 1.686 

-2.647 

asymptote 

3 

1.656 

-2.500 

TIME 

STEP 

1«6 

RISTAR 

a 2,570 

-1,797 

asymptote 

3 

2,489 

-1,675 

TIME 

STEP 

151 

RISTAR 

3 5,061 

-.664 

ASYMPTOTE 

3 

2,940 

-,595 

TIME 

STEP 

156 

RISTAR 

» 5,079 

,576 

ASYMPTOTE 

X 

2.944 

.575 

TIME 

STEP 

161 

Risr AR 

a 2.621 

1,730 

asymptote 

3 

2,500 

1,658 

TIME 

STEP 

1 6b 

RISTAR 

a 1.757 

2,617 

asymptote 

3 

1,675 

2,489 

TIME 

STEP 

171 

RISTAR 

a ,620 

5,101 

asymptote 

3 

,595 

2.940 

TIME 

STEP 

176 

RISTAR 

a -,6I5 

3,105 

asymptote 

3 

-.575 

2.944 

time 

STEP 

181 

RISTAR 

a -1,756 

2,628 

asymptote 

3 

-1.658 

2.500 

TIME 

STEP 

186 

RISTAR 

a -2,626 

1 , 746 

asymptote 

3 

-2,489 

1.675 

TIME 

STEP 

191 

RISTAR 

a -3,092 

,593 

ASY MPTOTE 

3 

-2.9U0 

.595 

TIME 

STEP 

196 

RISTAR 

* -3^082 

-.649 

asymptote 

a 

-2,944 

-.575 

TIME 

STEP 

201 

RISTAR 

3 -2,597 

-1,790 

asymptote 

B 

-2,500 

-1,656 



TIME 

INPUT r?(TlME) 







TIME 

INPUT 

3(TIME) 








0,000 .101 

.193 

,101 

-.000 

-.101 

-.193 

-.101 

0. 

1,000 

.707 

-,000 

-.707 

-1,000 

-.707 

-,’000 

,707 

l2bE-0^i 

-,lbb .ab8 

,827 

,702 

• 166 

-,968 

-.627 

-.702 

.I25E-09 

,929 

,929 

,383 

-.363 

-.929 

-.929 

-,363 

.383 

2!j0E-0« 

-.575 ,550 

1.355 

1,363 

.575 

-.550 

-1,353 

-1,363 

.250E-09 

,707 

1.000 

,707 

,000 

-.707 

•1.000 

-.707 

-.000 

J7SE-oa 

-1,IU3 ,329 

1,606 

1.995 

1.195 

-.329 

-1.608 

-1,995 

.375E-09 

.383 

,929 

,929 

.363 

-,5B5 

-.929 

-.929 

-.383 

500E-0a 

-1.750 -,159 

1.526 

2.516 

1,750 

• 159 

-1.526 

-2.316 

.500E-09 

,000 

,707 

1,000 

.707 

-,000 

-.707 

-1,000 

-,707 

b2bE-0<i 

-2.2b2 -,820 

1,102 

2,379 

2.262 

,620 

•1,102 

-2,379 

,625t-09 

-.585 

,383 

,929 

,929 

,383 

-.383 

-.929 

-.929 

7b0E-0<< 

-2,559 -1,525 

,397 

2,087 

2,559 

1.525 

-.397 

•2.087 

.750E-09 

-.707 

,000 

.707 

1,000 

,707 

-.000 

-.707 

-1,000 

87bE-04 

r2,537 -2.130 

-.975 

1.958 

2,557 

2.130 

.975 

-1,958 

.875E-09 

-,929 

-.383 

.383 

,929 

,929 

,363 

-,383 

-.929 

lOOE-03 

•2,171 -2,503 

-1,369 

,567 

2,171 

2.505 

1,369 

-.56 7 

.lOOE-03 

-1,000 

-,707 

,000 

.707 

1,000 

.707 

.000 

-.707 

UiE-03 

-1,979 -2,598 

-2,129 

-.956 

1.979 

2.598 

2,129 

,956 

•lliE-03 

-.929 

-,929 

-.383 

,363 

,929 

,929 

,383 

-.383 

12bE-03 

-.593 -2,2?9 

-2,602 

-1,956 

.593 

2,229 

2.602 

1.956 

•125E-03 

-.707 

-1,000 

-.707 

,000 

.707 

1,000 

,707 

.000 

138E-03 

.508 -1,552 

-2,703 

-2,271 

-.500 

1,55? 

2.703 

2,271 

.l3aE-03 

-,585 

-,929 

-.929 

-.303 

,363 

,929 

,929 

.303 

150E.03 

1,518 -.617 

-2.390 

-2.769 

-1,518 

.617 

2,390 

2.769 

,l50E-03 

-.000 

-.707 

-1.000 

-.707 

.000 

.707 

1,000 

,707 

163E-03 

2,329 ,998 

-1,695 

-2,895 

•2,329 

-,99R 

1,695 

2,895 

.163E-03 

,563 

-,303 

-.929 

-.929 

-.303 

.383 

,929 

.929 

175E-03 

2.807 1.980 

-.719 

-2.990 

-2,807 

-1,980 

.719 

2.M90 

,l75E-03 

,707 

-,000 

-.707 

-1,000 

-,707 

,000 

,707 

1,000 

10«E-Oi 

2,8b9 2,31b 

,906 

-1,791 

-2,869 

-2.316 

-.906 

1,791 

,166E-03 

,929 

,363 

-,383 

-.929 

-.929 

-.383 

.383 

,929 

200E-03 

2.993 2,818 

1,992 

-.708 

-2,995 

-2,818 

-1.992 

,708 

.200E-03 

1,000 

,707 

- , 000 

-.707 

-1,000 

-.707 

,000 

,707 

213E-03 

1,727 2,89b 

2.369 

.959 

-1,727 

-2,896 

•2,369 

-.959 

.213E-03 

,929 

,929 

,363 

-.383 

-.929 

-.929 

-,383 

,383 

22SE-03 

,b82 2.529 

2,699 

1.569 

-,682 

-2,529 

-2,899 

-1,569 

.225E-03 

,707 

1,000 

,707 

-.000 

-.707 

-1.000 

-.707 

,000 

238E.03 

-.989 1,763 

2.977 

2,998 

,989 

-1.7b3 

-2,977 

-2,996 

.238E-03 

.383 

.929 

.929 

.363 

-.383 

-.929 

-.929 

-.503 

250E-03 

-1.592 .711 

2,597 

2.962 

1,592 

-.711 

•2,597 

-2.962 

.250E-03 

.000 

,707 

1.000 

,707 

-,000 

-.707 

-1,000 

-.707 

263E-03 

-2,9b8 -,9b8 

1.607 

3,025 

2,968 

,966 

-1,807 

•3.023 

.263E-03 

-.383 

,383 

,929 

,929 

,563 

-.383 

-.929 

-.929 

27S6.03 

-2,972 -1,569 

.725 

2.M9 

2,972 

1.589 

-,7?5 

-2,619 

.275E-05 

-.707 

.000 

,707 

1,000 

.707 

-.000 

-.707 

-1.000 

288E-03 

-3.021 -2,978 

-,9B5 

1,795. 

3,021 

2.976 

.963 

-1.795 

.288E-03 

-,929 

-.363 

.383 

,929 

,929 

,335 

-.383 

-.929 

3U0E-03 

-2,603 -2,991 

-1,626 

,690 

2,603 

2.991 

1,626 

-.690 

.300E-03 

-1,000 

-.707 

,000 

,707 

1,000 

,707 

JOO 

-.707 

313E-03 

-1,778 -3.099 

-2.527 

-.529 

1,778 

3,099 

2,527 

,529 

.313E-03 

-.929 

-.929 

-.383 

,363 

,929 

,929 

,303 

-.303 

325E-05 

-.672 -2,629 

-3,090 

-1,675 

,672 

2.629 

3,090 

1,675 

.325E-03 

-.707 

•1.000 

-,707 

,000 

,707 

1,000 

,707 

-.000 

338E-03 

•59b -1,799 

-3,063 

-2,566 

-,596 

1 . 799 

3.083 

2,566 

,338E-03 

-.383 

-.929 

-.929 

-.383 

.383 

.929 

,929 

.583 

3S0E-03 

1,68b -,660 

-2.697 

-3,069 

-1,666 

,660 

2,697 

3,069 

.350E-03 

-,000 

-.707 

-1,000 

-.707 

,000 

.707 

1,000 

,707 

363E-03 

2.570 .597 

-1,797 

-3,068 

-2.570 

-.597 

1.797 

3.080 

.363E-03 

,383 

-.383 

-.929 

-.929 

-.383 

,383 

,929 

,929 

37SE-03 

3,061 1,695 

-.669 

-2.639 

-3,061 

-1,695 

,669 

2.639 

.375E-03 

.707 

-.000 

-.707 

-1,000 

-.707 

,000 

,707 

1,000 

388E^03 

3,079 2,589 

.576 

-1,770 

-3.079 

-2.589 

-.576 

1.770 

,388E-03 

,929 

,363 

-.363 

-.929 

-.929 

-,363 

,383 

,929 

auuE-03 

2,621 5,076 

1,730 

-,b30 

-2,621 

-3,076 

-1.730 

.630 

.900E-03 

1,000 

,707 

-,000 

-,707 

-1,000 

-.707 

.000 

,707 

«J3E-03 

1.757 5.093 

2,617 

,608 

-1.757 

-3.095 

-2,617 

-,608 

,913E-03 

,929 

,929 

,383 

-.383 

-.929 

-.929 

-.385 

.503 

a2bE-03 

,620 2,651 

3.101 

1,755 

-,620 

-2,631 

-3,101 

-1,755 

,925E-03 

,707 

1.000 

,707 

-.000 

-,707 

-1,000 

-.707 

,000 

,ai8E-03 

-.615 1,761 

3,105 

2,630 

,615 

-1.761 

-3,105 

-2,630 

.938E-03 

.583 

,929 

,929 

,303 

-.383 

-.929 

-,929 

-.363 

4S0E-03 

-1,756 ,617 

2,626 

3,100 

1,756 

-.617 

•2,628 

-3.100 

,950E-05 

,000 

,707 

1,000 

,707 

-.000 

-.707 

-1,000 

-.707 

ab3E-03 

-2,626 -.623 

1.796 

3,091 

2.626 

.623 

-1,796 

-5.091 

.963E-03 

-.385 

,383 

.929 

.929 

,505 

-.363 

-,929 

-.929 

M7bE-03 

-3.092 -1,767 

.593 

2,60b 

3,092 

1,767 

-,595 

-2.60b 

.975E-03 

-.707 

,000 

,707 

1.000 

,707 

-.000 

-.707 

•1,000 

i4b8E*03 

-3,002 -2,630 

-.699 

1,720 

3,062 

2,6 38 

,699 

-1,720 

.988E-03 

-,929 

-,383 

,383 

.929 

,929 

,383 

-,363 

-.929 

bOOE-03 

-2.597 -5,102 

-1.790 

,570 

2,597 

3.102 

1,790 

-.570 

.500E-05 

-1.000 

-,707 

.000 

,707 

1,000 

.707 

-,000 

-.707 



4 =- 

4 =- 


TIME 

G 9 EENS 

FUNCTION HR 



0.000 

25,065 

-,ooo 

-.000 

-,000 

0,000 

.IS7 

-1,661 

,078 

,280 

,517 

.520 

• 5H 

-2,505 

,163 

,370 

,398 

.000 

.071 

-2.650 

.332 

.589 

.395 

.397 

.6?8 

-2,606 

.508 

.301 

,316 

.310 

,705 

-2.008 

,d76 

.200 

,19a 

.185 


-1,999 

,800 

.in 

,0b0 

,050 

l.lOO 

-1,525 

,655 

-,01b 

-,ooo 

-.055 

1,257 

-1,061 

,815 

-.115 

1 08 

-.107 

X.tna 

-.751 

,666 

-.159 

-.111 

-,100 

1,571 

-,505 

,097 

-,10b 

-.069 

- , 058 

1.726 

-,595 

,200 

-.077 

-,007 

,003 

1,605 

-.567 

,076 

.030 

,051 

,050 

2.002 

-.579 

-,063 

.155 

,062 

.076 

2.199 

-.592 

-.160 

,2b7 

.077 

,060 

2.556 

-.580 

-.213 

.351 

,000 

.028 

2,513 

-.508 

-.197 

,590 

-.011 

•,01b 

2,670 

-.295 

- , 1 55 

,362 

-,05b 

-.007 

2,627 

-,25b 

-,110 

.332 

-.076 

-.050 

2,965 

-,191 

-,078 

.255 

-,ObO 

-.030 

5,102 

-.167 

, 067 

,167 

-,018 

.003 

3,299 

-.162 

-.075 

,085 

,052 

,039 

5,056 

-.170 

-.091 

,020 

,126 

.059 

3,613 

-,179 

-,106 

-.020 

,169 

,053 

3,770 

-.161 

1 1 1 

-.037 

.226 

,021 

5.927 

-.172 

-,103 

-,03b 

.237 

-.023 

't.080 

-.153 

-,067 

-,025 

,219 

-,Obl 

9.201 

-.152 

• ,067 

-.015 

,163 

-.075 

0,596 

•,110 

-.051 

-.006 

.159 

-.051 

0,555 

-.100 

-.000 

-,00b 

,097 

.011 

0,712 

-,10O 

-,005 

-.010 

,Obl 

,100 

0,869 

-.110 

-.052 

- ,02b 

,030 

.197 

5,027 

1 1 7 

-.060 

-.037 

.015 

,280 

5,160 

-,120 

• ,060 

-.003 

,002 

.529 

5,501 

-.116 

• , 060 

-.003 

-.007 

,53b 

5,098 

-.106 

-.051 

-,0 5b 

-.on 

,505 

5,655 

-.093 

-.000 

-.027 

-.008 

.202 

5,612 

-,085 

-.030 

-.017 

,000 

.170 

5,969 

-.078 

-.025 

-.012 

.020 

.100 

6,126 

-,079 

-.026 

-.012 

,050 

,05b 

6,265 

• , 065 

-.032 

-.018 

,077 

.031 

6,000 

-.091 

-.036 

-.028 

,099 

,025 

6.597 

-.090 

-.002 

-.037 

.Ml 

,030 

6,750 

-,092 

-.001 

-.002 

.no 

.037 

6,912 

-.080 

-.030 

-.039 

,095 

.039 

7,069 

-.075 

-.025 

-.027 

,072 

,030 

7,226 

-.067 

-.016 

-,006 

,006 

,022 

7,565 

-.060 

-.on 

,017 

,022 

,008 

7.500 

-.066 

-.013 

,001 

,00b 

•,000 

7.697 

-.072 

-,020 

.061 

-.001 

-.010 

7,650 

-.079 

-.029 

,075 

-.000 

-.009 


TIME 

GREENS 

Function hr 



e.on 

-.005 

-.036 

,075 

,004 

-.005 

6.168 

-.001 

-.037 

,069 

,009 

,000 

6,325 

-,070 

-.030 

,055 

.010 

,003 

6,062 

ObO 

-,Oia 

,03B 

,006 

,002 

6.639 

-.055 

,006 

,021 

,002 

-.003 

6, 79b 

-,052 

.027 

,009 

-,005 

-.009 

6.950 

-.057 

,095 

.001 

-.on 

-.014 

9,111 

- , 067 

,056 

-.001 

-,013 

-,01b 

9,268 

-.079 

,059 

,000 

-,012 

-.015 

9.025 

- , 0 87 

.055 

,003 

-.008 

-.011 

9,562 

-.065 

.046 

,005 

-.005 

-.008 

9.739 

-.071 

.035 

,006 

-.003 

-.005 

9,89b 

-,0OO 

,02a 

,004 

-,003 

-.005 

10,055 

-.009 

,015 

-.000 

-,00b 

-.008 

10,210 

,025 

,009 

-.005 

-.010 

-.012 

10.367 

.053 

.003 

-.006 

-.014 

-.015 

10,520 

,067 

-,000 

-.010 

-.015 

-.016 

10.661 

,066 

-.003 

-,010 

-.014 

-.015 

10,656 

,052 

-,005 

-.007 

-.010 

-.012 

10,996 

.029 

,000 

-.004 

-.007 

-.008 

11,155 

,005 

,007 

-.002 

-.005 

-.006 

11.510 

-.015 

,017 

-.001 

• , 005 

-.006 

U.0b7 

-.on 

• oai 

-,003 

-.020 

-.026 

11,620 

021 

,oaa 

-.007 

-.019 

-.024 

11.781 

0 2o 

,052 

-.on 

-.018 

-.020 

11,936 

- . 026 

.053 

-.015 

0 1 6 

-.017 

ir 095 

-.020 

,050 

-,015 

-.015 

-,01b 

12.252 

-.021 

.045 

-.012 

-.014 

-.016 

12.009 

-.020 

,038 

-.005 

-.014 

-.016 

12,566 

-.021 

,031 

.005 

-.015 

-.021 

12,725 

-.020 

.024 

,016 

-.017 

-.024 

12,661 

-.029 

,019 

.027 

-.019 

-.025 

15.058 

-.030 

,015 

.034 

-.020 

-.025 

15.195 

-.057 

.015 

,036 

-.021 

-.023 

15.352 

-.023 

.013 

,021 

-.020 

-.004 

13,509 

-.027 

,013 

.020 

-,016 

D05 

1 3,666 

-,030 

.013 

,017 

-.013 

-.005 

15,623 

-.052 

,015 

.013 

-.007 

-,OOH 

15.960 

-, 0 53 

.013 

,008 

.001 

-.012 

10.137 

-.052 

,012 

,005 

,000 

-.016 

10,290 

-.051 

.011 

-.001 

,015 

-.019 

10,051 

-.029 

,009 

-,005 

,020 

-.020 

10.606 

-.028 

.008 

-.006 

,022 

-,018 

19.765 

-.02/ 

.007 

-.009 

.022 

-.012 

10,925 

-.027 

,006 

-.010 

,019 

-.005 

15.060 

-.020 

,006 

-.009 

.014 

.007 

15.257 

-.050 

,006 

-.008 

,007 

,020 

15,590 

-.032 

.007 

-,007 

,001 

.031 

15.551 

-,03a 

.007 

-,006 

-.005 

,042 

15,706 

- .035 

.007 

-,00b 

-,009 

,049 



TIME 

GKEtiNS 

FUNCTION JR 

TIME 

GREENS 

FUNCUON JR 

0,000 

- 3,698 

- 1.520 

-.629 

8,011 

• 018 

-•129 

• •033 

• 157 

- 3,114 

- 1,248 

-,514 

8,166 

• 024 

••116 

-.035 

.310 

- 2.208 

-.851 

-,348 

8.325 

,023 

-.093 

-,029 

• ^71 

- 1,158 

-.420 

-.171 

8,462 

•on 

• •066 

-.023 

, 6 a 6 

-.078 

-.054 

-.024 

0,639 

-,013 

-,044 

•, P 17 

, 78 b 

,820 

,185 

.072 

6,796 

-.045 

-.032 

-.013 

,<»«2 

1,448 

,264 

,098 

0,954 

-.077 

-.029 

-.012 

1.100 

1.764 

,198 

,073 

9.111 

-.101 

•,032 

-.014 

1.257 

1,786 

,050 

.020 

9,266 

-.113 

•.037 

-.016 

1.414 

1,590 

-,108 

-,032 

9,425 

-.109 

•,039 

••016 

1.571 

1.267 

-.208 

-.059 

9.582 

-.093 

•,037 

-.016 

1,728 

,912 

-.203 

-.055 

9.739 

-,069 

-,030 

•,013 

1.885 

, b 01 

-,090 

-.024 

9,896 

-,046 

-,022 

-.010 

2.042 

,375 

,105 

,017 

10,053 

-.029 

-.016 

-,008 

2,199 

• 242 

.529 

,046 

10.210 

-.025 

•«014 

-.007 

2,356 

.184 

.525 

,050 

10,367 

-,027 

-,016 

-.007 

2,513 

.172 

,650 

• 024 

10,524 

-.037 

•,020 

-.009 

2.670 

.174 

.661 

-,020 

10,681 

-,046 

-.023 

••010 

2.827 

,170 

.625 

-.062 

10.838 

-.049 

-,024 

-.Oil 

2,985 

.152 

,510 

-.080 

10,996 

-.041 

-.020 

-.009 

3.142 

.123 

,373 

-.056 

11,153 

-.024 

-.015 

-. 00 ? 

3,?99 

,091 

,249 

.012 

11.310 

.000 

-.009 

-.005 

3,456 

.065 

.163 

• 114 

11,467 

,019 

-,015 

-.010 

3,613 

.051 

, 12 ! 

,?28 

11,624 

,037 

-.020 

••013 

3,770 

,046 

,113 

,325 

11,781 

.050 

-.023 

-•013 

3.927 

,052 

,122 

.363 

11,936 

.057 

-,023 

-•on 

4,084 

,055 

,131 

,390 

12.095 

• 056 

••020 

-.006 

4.241 

,054 

• 126 

,348 

12,252 

.046 

-,012 

-.001 

4,398 

• 046 

.111 

.273 

12,409 

,037 

-,002 

• 002 

4,555 

,037 

,065 

,188 

12,566 

.025 

•on 

,003 

4,712 

,027 

,059 

.117 

12,723 

.015 

,024 

.001 

4.869 

,020 

,043 

,075 

12,881 

• 008 

.036 

-.004 

5,027 

,019 

,039 

• 066 

13,038 

• 006 

,045 

-.010 

5,184 

.021 

,046 

,082 

13,195 

,007 

,050 

-.015 

5 j , 34 l 

,025 

,055 

.106 

13,352 

,004 

, .051 

-.010 

5,498 

,027 

.061 

.125 

13,509 

• 004 

,046 

• •003 

5,655 

,026 

.050 

.120 

13,666 

,005 

,042 

,006 

5,812 

• 021 

,045 

,089 

13,823 

,008 

,034 

.016 

5,969 

• 014 

.029 

,036 

15.980 

.011 

.026 

,025 

6.126 

,009 

.016 

-,027 

14,137 

• 013 

.019 

,033 

6,283 

.007 

• on 

-.066 

14.294 

.014 

• 014 

.059 

6,440 

,009 

,017 

-.127 

14,451 

• 013 

•on 

, 04 ? 

6,597 

• 013 

.031 

-.145 

14.608 

.011 

•on 

,042 

6,754 

.018 

.045 

-.139 

14,765 

.008 

.012 

,041 

6,912 

.019 

,051 

-.116 

14,923 

,005 

.014 

,039 

7,069 

,015 

.042 

-.086 

15.080 

,003 

,016 

.035 

7,226 

,009 

,016 

-.057 

15.237 

• 001 

,017 

.052 

7,363 

.005 

-,025 

-,057 

15,594 

• 001 

.016 

.028 

7.540 

-.001 

-.066 

-.028 

15,551 

.001 

,015 

,025 

7.697 

,001 

-.103 

-.027 

15,700 

.003 

,012 

• 022 

7,854 

.008 

-.125 

-,030 






1 


j:- 

vn 



Jfr 

ON 


TIME 

U.OOO 

.1 S 7 

GREENS 

25,465 

••187 

FU’^iCTION KR 
-.000 -.000 
.009 .03? 

-.000 

,035 

0,000 

.036 

TIME 

8,011 

GREENS 

-.009 

FUNCTION KR 
-.004 ,008 

• 000 

-.001 

. h \^ 

-.256 

,020 

.041 

,044 

,045 

8.166 

-.009 

-,0 C 4 

,008 

,001 

,000 

.«7 i 

-.292 

,037 

,043 

,044 

,044 

6.325 

-,008 

-.003 

.006 

.001 

• 000 

, b 2 Q 

-.294 

„056 

,038 

,035 

,055 

6.462 

-.007 

-.002 

,004 

.001 

• 000 

,78^ 

-, 260 

,075 

,027 

,022 

.021 

8 • 639 

-,00 b 

.001 

.002 

• 000 

-.000 

,9^2 

-.222 

,089 

,012 

,007 

,006 

6.796 

-,006 

,003 

,001 

-,001 

-,001 

1,100 

-.169 

,095 

-.002 

-.005 

- , 006 

6,954 

-.006 

.005 

.000 

-,001 

-,002 

l .^ S 7 

-.120 

.091 

-.013 

-.012 

-.012 

9.111 

-, 007 

,006 

-.000 

-,001 

-•002 

1.91 U 

-.081 

,076 

-.018 

-.012 

-.012 

9.266 

-.009 

,007 

,000 

-.001 

-.002 

l . i >7 l 

- ,056 

,055 

-,016 

-.006 

-.006 

9,425 

-.010 

.006 

,000 

-.001 

-•001 

1.728 

-,0> J 4 

.031 

-.009 

-.001 

,000 

9,582 

-.009 

.005 

,001 

-,001 

-,001 

1,865 

-.041 

,008 

,003 

, OOb 

• 006 

9,739 

-.006 

,004 

.001 

-.000 

-,001 

2, oa 2 

-.042 

-.009 

.017 

,009 

,008 

■9,896 

-,005 

,003 

.000 

-.000 

-.001 

^.199 

-.044 

-,020 

,030 

,009 

,007 

10,053 

-.001 

,002 

-.000 

-.001 

-•001 

2.356 

-.043 

-.024 

,039 

,004 

,003 

10,210 

.003 

,001 

-.001 

-.001 

-•001 

2.51 i 

•.039 

-.022 

,043 

-.001 

-.002 

10,367 

,006 

,000 

-.001 

-,002 

-,002 

2,670 

-,033 

-.017 

.042 

-.006 

-.005 

10,524 

,007 

-.000 

-,001 

-,002 

-.002 

2,627 

-.026 

-,012 

,037 

-.009 

• . 006 

10,681 

.007 

-,000 

-.001 

-.002 

-.002 

2,965 

• , 0? 1 

-.009 

,026 

-.007 

-.004 

10.636 

,00 b 

-.000 

-.001 

-.001 

-.001 

3,142 

-,019 

-.007 

.019 

-,002 

.000 

10,996 

.003 

,000 

-.000 

-.001 

-.001 

3.299 

-.018 

-.008 

.009 

,00 b 

,004 

11,153 

.001 

,001 

-.000 

-•001 

-.001 

3, uS 6 

-.019 

-,010 

,002 

,014 

,007 

ri .310 

-.002 

.002 

-.000 

-,00 l 

-,00 l 

3.613 

-.020 

-,012 

-.002 

,021 

,006 

11,467 

001 

,005 

-.000 

-.002 

-.003 

3.770 

-.020 

-.012 

-.004 

,025 

.002 

1 1,624 

-.002 

,005 

-.001 

-.002 

-.003 

3,927 

-,019 

-,0 i 1 

-,00 U 

,026 

-.003 

11.781 

-.005 

,006 

-.001 

-,002 

-,002 

9,064 

• , 0 1 7 

-.010 

-.003 

,024 

-.007 

11,938 

-.003 

,006 

-.002 

-,002 

-.002 

9,241 

• , 0 15 

-.007 

-.001 

,020 

-.006 

12.095 

-.003 

,006 

-.002 

-.002 

-.002 

4,398 

-.013 

-.006 

-.001 

,015 

-.006 

12,252 

-.002 

,005 

-.001 

-.002 

-.002 

4,555 

-,012 

-.005 

-.001 

,011 

,00] 

12,409 

-.002 

.004 

-.001 

-,002 

-.002 

4.712 

-.012 

-,005 

-.002 

,007 

,cn 

12,566 

-,002 

,003 

• 001 

-,002 

• ,002 

4.869 

-,012 

-.006 

-.003 

,004 

,022 

12.723 

-,003 

.003 

• 002 

-.002 

- t 003 

5,027 

-.013 

-.007 

-,004 

,002 

.031 

12,881 

-,003 

,002 

• 003 

-.002 

-•003 

5,164 

-.013 

- , 007 

-.005 

,000 

.037 

13.036 

-.004 

,002 

.004 

-. OO ? 

-.003 

5.341 

• , 0 1 3 

-.007 

-.005 

-.001 

,037 

13.195 

-.004 

,001 

• 004 

-.002 

-,003 

5,498 

-.012 

-,006 

-.004 

-.001 

,034 

13.352 

-,003 

.001 

• 002 

-.002 

-,000 

5,655 

-, 0 1 0 

»«,004 

-,003 

-.001 

.027 

15,509 

-,003 

,001 

• 002 

-•002 

•tOOO 

5,612 

• ,009 

-,003 

-.002 

,000 

,019 

13, 666 

-.003 

,001 

• 002 

-.001 

-.001 

5,969 

-,009 

-.003 

-.001 

,003 

,012 

13,623 

-.004 

.001 

.001 

-.001 

-.001 

6,126 

-,009 

-,003 

-.001 

,00 b 

,006 

13,960 

-,004 

,001 

• 001 

,000 

-.001 

6,263 

-,009 

-,004 

- , 002 

,009 

.003 

14,137 

-.004 

,001 

.000 

,001 

-.002 

6,440 

-.010 

-,004 

-,003 

,011 

.003 

14,294 

-.003 

,001 

-•000 

,002 

-.002 

6,597 

-,010 

-,005 

-.004 

,012 

,003 

14.451 

-.003 

,001 

-•001 

,002 

-.002 

6,754 

-.010 

-.005 

-,005 

,012 

.004 

14,608 

-,003 

,001 

-•001 

,002 

-.002 

6,912 

-.009 

-.004 

-.004 

,011 

.004 

14.765 

-.003 

,001 

-.001 

,002 

-.001 

7,069 

-.008 

- , 003 

-,003 

,008 

,004 

14,923 

-.003 

,001 

-•001 

• 002 

••000 

7,226 

-,007 

-,002 

-.001 

,005 

,002 

15,060 

-.003 

.001 

-.001 

,002 

.001 

7.363 

- , 007 

-.001 

,002 

,002 

,001 

15,237 

-.003 

,001 

-•001 

.Ot)! 

,002 

7.540 

-,047 

-,001 

,005 

,001 

-.000 

15, S 94 

-.004 

,001 

-.001 

.000 

.003 

7.697 

-,006 

-,002 

.007 

-.000 

-,ooi 

15,551 

-.004 

,001 

-.001 

-.001 

,005 

7.654 

-.009 

-,003 

,008 

-,000 

-.001 

15,706 

-.004 

,001 

-•001 

-.001 

,005 


TIME 

PRESSURE PCTIN^E) 

FOR DISCRETE VALUES 

OF Y 




0. 

■ 526E+ 06 

.372E+C6 

-,7b0E-09 

-.372E+06 

- . 526E*0b 

-•372E406 

-.B22E-03 

.372E+06 

,l25e-0tt 

.S«2E*06 

.S2St+06 

,20SE+06 

-.238E4C6 

-.5y2E406 

-.S26E406 

-,20SEt06 

,23BEt0b 

.2S0fc-0a 

•620E+06 

•701E*06 

,U26E*0b 

-.155E+06 

-,b20E+0b 

-,7atE406 

•.tZRE+Ob 

.135E+06 

.375E-0U 

.uuOE+06 

, 7bbE*06 

.6«aE+0b 

, I y«E+06 

*,au0E+0b 

-.766E406 

-.6U«E406 

•,l«yE406 

,i>00E«0tt 

,R2oE+C5 

.599E406 

,7SSE+0b 

,U6ttE+06 

-,926L+0b 

-,S99E+0b 

-,755E+0b 

..«6eE+06 

•O25E-00 

-,577Et06 

.235£*0b 

,709£+0b 

, 7b7C + 0b 

,377h+0b 

-,235t+06 

-,709£+0b 

• , 7b7E*06 

, 7boE-0*i 

•.fl5RE+0b 

•• 256E + 06 

,a7bt+0b 

.930E+06 

,639E+0b 

,2SbE+0b 

•,«7bE*0b 

-,930E+Ob 

,.i7^£-0J 

-,121£*C7 

-,798r*c6 

.798E+05 

.911£ + (j6 

.121E+07 

,79et+0b 

-,79PE+05 

-,9llt+Qb 

.lOOc-OS 

-.l>SE+07 

12bt*07 

•, u27E + 0b 

,673E+0b 

.158E+07 

.123E+07 

.U27E+Cb 

673E+06 

,lliE-03 

-.131E+07 

lb0E*07 

-,9S5Et06 

.2b2E+06 

.131E+C7 

. 160E407 

,9b5tf06 

« , 252E + 06 

.125E-03 

Rd3E* 06 

-.:69E407 

lulc 407 

-,30lEtOb 

,983t*06 

,lfa9E+07 

,U1E*07 

.301E+06 


-.u‘j7E + 06 

-.151E+07 

-,169£+07 

•, ridOE + 06 

,yu7E+06 

,1S1E407 

.169E+07 

.8e0E+06 

. I^>0h*03 

•222E+06 

-,107t+07 

1 7«E+07 

-,139E+07 

-,222t+06 

.107E+07 

.17UE407 

.139E+07 

• 1 &5E-0i 

,R0y£+0o 

-,u35£ +06 

•a^2E + 07 

-.’7U + 07 

•,900E+Ob 

.U35E+06 

,l52Et07 

.17IE+07 

.l75t-05 

.lyRE+07 

.312E+06 

lOyE + 07 

-,1 79E+07 

-.1U9E+07 

-,312E+0b 

.lOuE+07 

.179E+07 

.iebE-05 

. 186E + 07 

, lOUh+07 

-,3buE+0b 

•,156Et07 

-,16bE+07 

109E + 07 

,38uE*0b 

.lb8E+07 

,2o0t-03 

. lRyit07 

. Ib3t*07 

,3b9t+06 

-,niE + 07 

-,19UE^07 

•.163E+07 

369E + 06 

.UlE + 07 

.21fL-0i 

,:7?£*0/ 

,l99tT.7 

, 1 OOE + 07 

-,y39t+0b 

172E + 07 

..199E+07 

- . 1 0 9 E + 0 7 

.«39E*06 

,22bt.-03 

.120E+07 

.203e407 

, 167E+07 

,336E+06 

120E + 07 

-.203E407 

-.tb7£407 

-.336E+06 

,23dt-0i 

, y ?0E + 06 

♦17btf07 

,201E+07 

.109E*07 

-,u70E+0b 

-.175E+07 

-.201E407 

1096407 

,2b0t-03 

•• 3S8E + 0b 

. 1 ieE*07 

.203E+07 

.169E+07 

.358E+06 

-,118Et07 

-.203E+07 

-.169E+07 

,2e5E-05 

liyE*07 

,yi9E+0b 

.I73tf07 

.203E407 

,naEt07 

-,ill9E + 06 

•,173fc*07 

-,203fct07 

,27b£-05 

-.176E+07 

-,«2St+0b 

,1 lbE + 07 

.206E+07 

,17fat. + 07 

,y2bE+0b 

• , n bE*07 

-.206E+07 

,2BdE-0i 

-,2llEi07 

-.122E+07 

,391Et06 

, 177E + 07 

.2nE + n7 

, J22E + 07 

-.391E+06 

-.177E+07 

,300£-05 

-.2l3Et07 

185E + 07 

-.y53E+06 

.U9E + 07 

,213E»07 

,lB3t+07 

,y53E+06 

-.M9E407 

.3:i£-03 

-,132t+07 

-,2lbt407 

129E + 07 

,ul2E+06 

,182Et07 

.216E+07 

.12y£+07 

-,ul2t+06 

•32SF-03 

-,122E*C7 

•.23 7M07 

1 8yE + u7 

-,yul t+06 

.122F+07 

.217E+07 

,ieyE+07 

,uy lE+06 

, iist- 0 i 

-,Ul<?E + 06 

l83t4Q7 

-.217E+07 

-,12yE+07 

,yi9E+06 

.103E+O7 

.217E+07 

.12^6407 

.3b0E-03 

,y56E*06 

lPlt+07 

-.216E+07 

-,l85E+07 

-,«56E+0b 

.121E+07 

.216E+07 

.18SE+07 

, 3b3t-03 

,127c*07 

-,392t40b 

-.1 «2Ef07 

-.218E+07 

-.127E+07 

.392E+06 

,lfl2fc+07 

.2IHE407 

.i/bt-Oj 

, lft'?Et07 

, u 9 0 E + 0 6 

-.1 19E + 07 

-.21 8E+07 

-.189E+07 

-.U9UE+06 

.M9E + 07 

.210E4O7 

•3SbE-03 

,222E+07 

.130F*07 

-,378Et0fa 

-,18«E+07 

-.222E+07 

-.130E+07 

.376E+0b 

,1396+07 

•UO^E-03 

.221E+07 

.191E+07 

.501E+06 

-.121E+07 

-.221E+07 

-.191E+07 

-,50U*06 

.121E+07 

.‘iiiE-OS 

. !dScfC7 

.223E407 

,l3!E+07 

-,383E*06 

-.lPSE+07 

-.223E+07 

-.131E+07 

•383E+06 

, £ b E - 0 3 

,120£*07 

,221 £*07 

,l9?E+07 

,b03E+Ob 

-.120fc+07 

-.221E+07 

19^Ef07 

50 3E + 06 

*‘4ib£-03 

.37CL*Ob 

, 1 y E + 0 7 

•223E+0? 

, 131E + 07 

•.370E+06 

-.1 5yfc+07 

-,225£+07 

-.1316+07 

.ii50t-03 

• ,5?y E ♦Ob 

, 1 19E + 07 

,220Et07 

.193E+07 

,52'JE + Ob 

-.119E+07 

-.220E+07 

-.1936+07 

,iio3t-03 

l3yE^07 

• 5 y 7 1 ♦ 0 6 

,lH3L+07 

,22yE+07 

.13yt*07 

-,3a7E*06 

-.163E*07 

-.229E+07 

,« /b£-Oi 

-.195E-J7 

•,5ybE+06 

,M8E+07 

.221E+07 

.195E+07 

,5ybE+06 

-,ll8Ef07 

-.221E+07 

.u8at-03 

-,22bEt07 

1 3EE + 07 

,3*40Et0b 

,lfiyE+07 

.226E+07 

.136E+07 

•, 390E t Ob 

-.1896+07 

•500t-03 

-,222E*07 

-,i96E*07 

-.5S0E+06 

, 1 18E + 07 

.222E+07 

. 196E + 07 

.SbOEf-Ob 

-.1186+07 


4 =- 


Jr- 

00 


TIME 

upstream component of pressure 





0. 

.S00Ef06 

*3b«E+06 

-.65SE-09 

-.354E+06 

-.500E+06 

-.3546+06 

-.8Q9E-08 

.354E+06 

•l25E-oa 

.«3aE*06 

,U39E+06 

.187E+06 

-.175E+06 

-.434E+06 

-.43RE+06 

-.187E+06 

,1756+06 

,2bOE'-Ott 

.330E+06 

,«7?E+06 

.338E+06 

.595E+0U 

-.330E+06 

-.4726+06 

-.338E+06 

-.5956+04 

,37bE»oa 

, 169E + 06 

.«29e+06 

.«3eE+06 

.190E406 

-,l69C+06 

-.429E+06 

-.438E+06 

-,1906*06 

,500E-0tt 

-,l«SE+05 

.321E+06 

,«69c+06 

.3426406 

,1456+05 

-.321E+06 

-.4696+06 

-.342E+06 


-.199E+06 

.IblE^'Ob 

,u27t+Ob 

,4tt2t40b 

,1996+06 

-.1616+06 

-.4276+06 

-.4426+06 

,7b0E-0a 

-,3blE+06 

-.2U0E+05 

.317E+06 

,4736406 

,3516+06 

,2406+05 

-.3176+06 

-.475E+06 

,a/5E-0« 

-.«5lE+a6 

-,207t+06 

. 158E+06 

.431E406 

,4516+06 

.207E+C6 

-.1S8E+06 

-.4316+06 

.IQOE-Oi 

•.ufiOE^Ob 

-,359E+0b 

-.273E+05 

.320L406 

,4806+06 

,3596+06 

.2736+05 

-.3206+06 


-,i<37E + 06 

-,4S7E+0b 

-.210C+06 

. 161E406 

.437E+06 

,4576+06 

.2106+06 

-, 1616 + 06 

.l25E-0i 

-,52iiE + 06 

-,*jbbE + 06 

•,3blE+06 

.258E40S 

,3246+06 

,4656+06 

,3616+06 

.2586+05 

.libt-OS 

-,163E^06 

-,y39t40fe 

-,a56E+06 

-.2096406 

.165E+06 

,4396+06 

.458E+06 

.2096+06 

,lbOt-Oi 

,?blEt05 

-.32bfc+0b 

•.apSEtOb 

-.361E406 

-.251E+05 

.3256+06 

.485E+06 

.5616*06 

,lo3E-0J 

,209cf06 

-•lb?t+06 

-,a3‘5fc + 06 

-,4b8E406 

•.209E+06 

,1626+06 

,4596+06 

.4586+06 

•175C-03 

, 562E + C6 

,2bOf +05 

-.32SE+06 

-.4866+06 

-,3626+06 

2606^05 

.3256 + 06, 

, a8bE*06 

,iaa£-03 

,ub9£f06 

,210E*0b 

• . 1 62F +06 

-.4396 +06 

-.4596+06 

-.2106+06 

, 162E + 06 

,4396+06 

,2aot.o3 

^ii67L*0b 

.353Et06 

.265E+ns 

3266+06 

-.4876+06 

3636 + 06 

-.265E+C5 

.3266*06 

.213E-03 

.a^oE+Ob 

,«JblE + 06 

.211E+06 

-,1626+Ob 

- .440F + 06 

-.4616+06 

-.21 lE+Ob 

.1626+06 

,225E-03 

.52b£+06 

,a8BE+C6 

.369E+06 

,2b4E+05 

-.326E+06 

-.4886+06 

-.364E+G6 

-.2646+05 

.23CE-03 

.162E+06 

.UuOE+06 

,46lE+ab 

.2116+06 

-.162E+06 

-.4406+06 

-.4616+06 

-.21 IE+06 

,250£-03 

-.272E+05 

.325E+06 

.487E+06 

.364E+06 

.2726+05 

-.3256+06 

-.4876+06 

-.3646+06 

.263E-03 

-,2t2E+06 

.Iblt+Ob 

,4UOE+Ob 

.461E+06 

,2126+06 

1616 + Ob 

-.440E+06 

-.4616+06 

•275E-03 

3buE+06 

-,272E+05 

,325E+0b 

,4876+06 

,3646+06 

,2726+05 

-.325E+06 

-.4876+06 

•2tf0t-03 

-,4olE*C6 

-.212E+06 

, lblE + 06 

,4406+06 

,4bl 6+06 

,2126+06 

-.1616+06 

04CE + 06 

, 5 O 0 E .-05 

-,^67E*06 

3buE + C6 

-.272E+05 

.325F+C6 

.4876+06 

.36UE+06 

,2726+05 

-.3256+06 

,5iiE-0i 

-.u^OE*Ob 

-.ublE+Ob 

-.212E+06 

.1616*06 

,4406+06 

,4616+06 

,2126+06 

-.1616+06 

.325E-03 

-.325E+06 

-.y87E+0b 

-.56UE+06 

2726 + 05 

.3256+06 

.4876+06 

,3646+06 

.2726+05 

, 33t3c*0 3 

-,lbU + Ob 

•,tiiiOE + Ob 

-.461E+06 

-.2126+06 

,1616+Ob 

.4406+06 

.4616*06 

.2126+06 

.35OE-03 

•272Ef0b 

-.52SF+06 

-.487E+0b 

-.3646+06 

-.272E+05 

.3256+06 

.4076+06 

,3646+06 

.3636-03 

•212E+06 

-.lblE+06 

-,4a0E+06 

-.4616+06 

-.212E+06 

.1616+06 

, UUOE + 06 

, 4ol E + 06 

.37-DE-03 

,3b«E+06 

.272E+OS 

-.325E+06 

«.«87F+06 

-,5b4t+06 

2726 + 05 

,3256+06 

,4876+06 

.38ct-03 

.ut) I £ + 06 

.212E+06 

1 blEfOb 

-.440E+06 

-.4616+06 

-.2126+06 

.1616+06 

,4406+06 

.aoOE-03 

,aS7t *06 

,3bot+Gb 

.272E+05 

-.3256+06 

-.4876+06 

-.3646+06 

-.272E+05 

.325E+06 

.Jl 3E-0S 

,uaOE+06 

,ab\E+Ob 

.212E+06 

-.lblE+06 

-.4406+06 

-.4616+06 

-.2126+06 

.1616+06 

.uasE-os 

,325l+06 

,ab7t+C6 

.3o<iE+06 

,2726+05 

-.3256+06 

-.4676+06/ 

-.5646+06 

-.2726+05 

.u3ae-03 

, Ibl E + 06 

,<i^i0h + 06 

,4olE+06 

.2126+06 

-.1616+06 

-.4406+06 

4bl E + 06 

..212E+06 

.u50fc-03 

-,272E+0S 

.3256+06 

.467E+06 

,3646+06 

.2726+05 

-.3256+06 

-.4076*06 

-.3646+06 

.4o3t-03 

-.212E+06 

.161 t+06 

,a40E+06 

, 4b 1 E+06 

,2126+06 

-.1616+06 

-.4406+06 

-.461 e+06 

.•*7b£.i)3 

364E+06 

-.272f +05 

.325E+06 

.4P7E+06 

.364E+06 

,2726+05 

-.325E+06 

-,457t +06 

,ttttttE-03 

-.ttblEiOb 

-.212E+06 

,16lE+0b 

.4406+06 

, 4616 + 06 

.2126+06 

-.161E+C6 

-.4406+06 

.500E-03 

-,a67t+06 

-,3biiE + 06 

-.272E+05 

,3256+06 

,4676+06 

,3646+06 

,2726+05 

-,325t+0b 
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